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Journal 
The Franklin Institute 


Devoted to Science and the Mechanic Arts 


Vol. 272 OCTOBER, 1961 No. 4 


ON NONLINEAR OSCILLATIONS IN ELECTROMECHANICAL 
SYSTEMS * 


BY 


Y. 


ABSTRACT 


Part I of the paper gives the representation of an electromechanical system by 
the dynamic equation of Lagrange in the holonomic reference frame. The equation 
of voltage and the equation of torque are then obtained from Lagrange’s equation. 

Part II is concerned with the transformation of the system from the holonomic 
frame to the nonholonomic frame. Lagrange’s equation is generalized to the Boltz- 
mann-Hamel equation. Theory of hunting is presented by using a uniformly rotating 


reference frame. 

Part III deals with the transformation of the system from the holonomic frame to 
the quasi-holonomic reference frame in which the d, g components (Park variables) 
and the f, components (Ku variables) are developed. From the equation of voltage, 
the impedance matrix is obtained by considering w, the angular velocity, as a con- 
stant. The equation of torque is, however, nonlinear. If w is considered as the 
derivative of the angular displacement 6, the “fourth dimension" in an electro- 
mechanical system, the system of equations are nonlinear. A simple four-dimensional 
formulation thus represents the dynamic characteristic of an electromechanical 
system. Four simultaneous first-order nonlinear differential equations are explicitly 


given. 
Equations of hunting in the nonholonomic frame and equations of voltage for 
transient analysis in the quasi-holonomic frame, are given in Appendices I and II, 


respectively. 


* Presented at the International Symposium on Nonlinear Vibrations, Kiev, USSR, 
September 12-18, 1961. This symposium was organized under the auspices of the Inter- 
national Union of Theoretical and Applied Mechanics, of which the author has been an elected 
member of the General Assembly since 1946. 

1 The Moore School of Electrical Engineering, University of Pennsylvania, Philadelphia, Pa. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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Part I 
1. THE EQUATION OF LAGRANGE 


An electromechanical system can be represented by the dynamic 
equation of Lagrange: 


d f{ aT ol oF 
where JT denotes the total kinetic energy function and F denotes the total 
energy loss in the electromechanical system, x7 denotes the coordinates, 
7 their derivatives with respect to time, and e, denotes the electro- 
motive forces or torques. The total kinetic energy function is given 
by the quadratic form 


T = (2) 


which comprises both the energy stored in the magnetic field W,, and 
the energy stored in a rotating body W., where 


Wa = Tm = (3) 


W. = = (4) 


The total energy loss (the Rayleigh dissipation function) in the electro- 
mechanical system comprises both the electrical loss F, and the friction 
loss in the rotating mechanical system F;. In quadratic form these are: 


(5) 


1 
3% 


= 


(7) 


We see that the symbol i*, denoting the derivative of x* with respect to 
time, represents either the electric current 7” or the angular velocity w". 
Hence the symbol x*, denoting the coordinate, represents either the 
electric charge g” or the angular displacement 6". Note that x7, where 
the index y includes all values of a and 8, is a contravariant tensor, and 
e, is a covariant tensor. The indices (a, 8) in Eqs. 2 and 5 include the 
electrical coordinates (m,m) and the mechanical coordinates (u, v). 
Since we have introduced two sets of indices (m, m) and (u, v), there is 
no need to use the symbols g, 8, 7, w, L, J, R, and B. Thus 


) 
ann = mn | 


= 


Quy 
r (8) 
Tan n 
| 
luv = 


3 
| 
[J. 
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| 
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Denoting «7 by 77, Eq. 1 can be rewritten as 


d oT oF 


We note that e,, denotes the electromotive forces and e, denotes the 
torques. 

The coordinates, whether they are electric charges x"(g") or me- 
chanical angles x“(@"), are rigidly attached to the conductors or the 
frame. This reference frame is termed holonomic, and Lagrange’s equa- 
tion is valid in a holonomic reference frame. 

In a rotating machine (electric generator or motor), there are con- 
ductors on the stator as well as on the rotor. So 7™ would include both 
stator and rotor currents, and e,, would include both stator and rotor 
voltages. We further note that the inductance dn, = Lm, would in- 
clude self and mutual inductances. If the effect of saturation is taken 
into account, the inductance would be a function of the current 7”. 
The mutual inductance between the stator and the rotor is a function 
of the angle x“ or 6". In a salient-pole machine, there is also second- 
harmonic variation in the mutual inductance between the stator phases. 
Thus ay, = Lm, is a function of the coordinate x“. 

Substituting Eqs. 2 and 5 into Eq. 1 or 10 gives 


dt ax® 2 


Now if dag is not a function of 7 (neglecting saturation), the fourth and 
fifth terms are dropped and Eq. 11 simplifies to 


where the Christoffel symbol is defined by 


= 2 \ axe axe J (13) 


: 
) 
2” 
x on 6 u 
= 
ky 
3 
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Equation 12 was first given by Kron (1)? for the electromechanical 
system and was derived in Bewley’s recent book (2). 


2. THE EQUATION OF VOLTAGE 

If we limit e, to e,,, which denotes the electromotive forces or vol- 

tages, and assume that a,,, is a function of x" (or @) only, the equation 
of voltage is given by 


. di” 
+ Amn + = (14) 
dt ax" 


Although in general x“ would denote 6, 62, ---, in practice we may need 
only to consider one angular displacement @ and its corresponding angu- 
lar velocity w. Thus Eq. 14 simplifies to 
dt" | mn . 
Ca = + Gan + 36 wt". (15) 
In order to get a complete solution for the performance of the electro- 
mechanical system, the equation of voltage has to be solved simul- 
taneously with the equation of torque. In other words, the complete 
system represented by Eq. 11 or 12 has to be analyzed as a whole. 
Even in the non-saturation case, Eq. 15 is nonlinear in the sense 
that the coefficient for 7" in the third term of the right-hand member 
involves another dependent variable 7. = 6 = w. Unless we assume 
that w is a constant we cannot solve for a set of currents which satisfies 
Eq. 15. In the transient analysis of rotating machinery, we often do 
assume that during the short interval of electrical oscillations, w re- 
mains constant, and hence Eq. 15 reduces to a system of varying- 
parameter equations. However, it is well known that the system 
cannot be analyzed in a holonomic reference frame. In a synchronous 
machine study, the stator or armature currents (and voltages) have to 
be referred to the rotating reference frame (3, 4, 5). In the transient 
study of induction motors, one can choose to refer the stator quantities 
to the rotor or the rotor quantities to the stator. 
Maxwell's equation of voltage (Faraday’s law) as applied to a 
holonomic reference frame is given as 


d 
Rinnt (Linnt") (16) 


Cm = n 1 + 


where ¢ = L,»,t" denotes the flux-linkages. Let La, be both a function 
of the current 7” and the mechanical angle 7« = 6". Then 


dt ) dt + 1 ( ot 4 (17) 


? The boldface numbers in parentheses refer to the references appended to this paper. 
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Substituting Eq. 17 into Eq. 16 gives 
. dL mn ( ) 


(18) 


= + Las + + dt 


dt 


This expression reduces to Eq. 15 for neglecting saturation and limiting 
6“ to one angular displacement 6. The expression (dL,,/90) is some- 
times termed the motion matrix. 


3. THE EQUATION OF TORQUE 


If we limit e, to e, (or r), which denotes the torque, and assume that 
dm, is a function of x" (or @) only, the equation of torque is given by 


dv 1 Od inn 


+ >, 
d 


te (19) 


Since we are considering only one machine, there is only one torque 
and Eq. 19 simplifies to 


(20) 


The superscript of @" can be dropped if there is only one angular dis- 
placement involved. 

The power delivered to a motor or developed from a generator is 
given by the general expression 


P = ént”. (21) 


Substituting Eq. 15 into Eq. 21 gives 


On the righthand side of Eq. 22, the first term denotes the ohmic loss, 
the second term represents the power associated with the inductance 
voltage drop d,,.(di"/dt), and the third term denotes the “rotational 
power” which is associated with (w). In the case of motor, the rota- 
tional power represents the useful output, while P = e,,7" represents 
the input. Dividing rotational power by » gives 


mn 
T= (23) 


Since dn = Anni" according to Eq. 16, Eq. 23 can be rewritten as 


(24) 


257 
: 
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d 
Bo+ J = 7, 
u 
dt 2 006 
j 
: 
Ohm . 
= — 4m 
06 
: 
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This torque expression was derived by the author (5, 6) for a three- 
phase synchronous machine as 


(25) 


where 74, 72, and 7° denote the phase currents of the armature (stator), 
and $4, ¢s, and ¢c denote the phase flux-linkages. As we have noted, 
all the quantities in Eqs. 23, 24, and 25 refer to the holonomic reference 
frame. These expressions agree in form with the third term on the 
left side of Eq. 19 or 20. As a,,, is a function of 6 in the case of a general 
rotating machine with salient-poles, it is easily seen that the expression 
for the electromagnetic torque given in Eq. 19 gives the same result as 
in Eq. 25 with a proper sign change. 


Part II 
4. THE NONHOLONOMIC FRAME 


From both the theoretical and the practical points of view one 
should be able to transform from a holonomic system to a nonholonomic 
system in which the reference axes are no longer fixed to the conductors 
(on the stator and on the rotor) or the frame. A nonholonomic refer- 
ence frame may be termed a free frame. We have mentioned in Section 
2 that in the transient study of induction machines, one can choose to 
refer the stator quantities to the rotor or vice versa. However, if one 
is interested in connecting an induction machine with a synchronous 
machine, the rotor of the induction machine will not run at the syn- 
chronous speed which is the speed of the synchronous machine. As 
shown in (3, 4, 5, 6), both the stator and rotor quantities have to be 
referred to a synchronously rotating reference frame before interconnec- 
tion with a synchronous machine can be made. Physically, the arma- 
ture terminals of a synchronous generator and the stator terminals of 
an induction motor can be easily connected together through a constant- 
parameter network or a distributed-parameter transmission line. In 
fact one can make the connection in the stationary frame and transform 
these stator (or stationary frame) quantities to a synchronous rotating 
reference frame. Taking the induction machine as a simple electro- 
mechanical system we can see that it is desirable to refer both the 
stator quantities and the rotor quantities to a synchronous rotating 
reference frame. This synchronous rotating frame is thus a non- 
holonomic frame. 

Let the transformation tensor between the old and new coordinates 


be Then 
(26) 


: 
A B Od Cc 
06 06 06 oe 
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where 7?’ denotes the new currents or velocities. Since the power or 
energy in the system is invariant regardless of the particular set of 
variables, we get 
P = = (27) 
It follows that 
ey = Cy ey (28) 


where e,, denotes the new voltages or torques. When C,* is non- 
singular, it is possible to find its znverse. In general, however, C.* may 
be singular, and the invariancy of power calls for a relation given in 
Eq. 28 in which C,,7 is to be interpreted as the transpose of C..* given 
in Eq. 26. 

Transforming Eq. 12 gives 


= Cye, = (Cy 
di’ 


(<= iv’) all 


+ Cy[oB, (29) 


+ (Cy 


Thus Lagrange’s equation can be generalized to the Boltzmann-Hamel 
equation 


_d 
er at \ ax’ 


where 
T’ = (31) 
FY = (32) 
= (33) 


= Cp? 
In Eqs. 33 and 34, C.-* is the same transformation tensor as that given 
in Eq. 26, and C36 corresponds to C,,7 in Eq. 28. In matrix notation, 
we have in general 

Z' =C.ZC (35) 
where Z denotes the old (original) impedance matrix, and Z’ denotes 


the new impedance matrix. Substituting Eqs. 33 and 34 into Eq. 29 
and introducing the affine connection give 


die’ 
= Qy'a’ a + 48 (36) 


: 
: 
? 
: 
2 
: 
: 
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where the —° -* connection is obtained by identifying the sum of the 
third and iourth terms on the right side of Eq. 29 with the third term 
on the right side of Eq. 36, or 


Cy? + [aB, (37) 


4! = Aya 


Note also that 


acy’ 
ax ax? 


II 


(38) 


( 


5. THEORY OF HUNTING 


Kron (7) first suggested the use of a uniformly rotating reference 
frame for the study of the hunting of one or two machines. An inde- 
pendent derivation of Kron’s result was obtained by the author in the 
summer of 1951 while working as a consultant for the General Electric 
Company. The author submitted a comprehensive paper (6) including 
the hunting equation to the Institution of Electrical Engineers, London, 
in February 1952. At the urge of the author, Kron presented his ‘‘New 
Theory of Hunting’’ before the A.I.E.E. Summer General Meeting, 
June 1952. This paper later appeared in the A.J.E.E. Transactions (8). 
To quote from Kron: 


If, however, the hunting equations of a synchronous machine are estab- 
lished in such a manner that the armature reference frame merely rotates 
uniformly with respect to the stationary armature and does not oscillate, 
then the equations and equivalent circuit of the series-connected trans- 
mission lines become comparatively simple, unencumbered by additional 
generators in each mesh. . . . Any number of synchronous and induction 
machines and transmission lines then may be connected together as if they 
were merely stationary impedances. 


It may be pointed out that in establishing an equivalent circuit it is 
necessary to use the forward and backward components or Ku variables 
(3, 5, 9) in the uniformly rotating reference frame instead of the direct 
and quadrature components, or Park variables (10). 

In order to simplify the discussion, we start from Eq. 16 for the 
holonomic frame with p = d/dt. Thus Eq. 16 becomes 


en = Rani" + Pbm. (39) 


Transforming Eq. 39 to the nonholonomic frame (uniformly rotating 
frame) gives 


Cm? = Rant” + (40) 


. 
| 
a 
2 
: 
A 
: 
: 
‘ 
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where p’ denotes the absolute derivative and can be given by 
pb’ = p + (41) 


This symbol y denotes the “rotation”? tensor which comprises the 
elements +7 and —j in connection with the f (forward) and 6 (back- 
ward) components. When only one machine is hunting, (6) repre- 
sents the actual angular velocity of the rotor, and (p6.) = w represents 
the fundamental angular velocity (frequency) of the generating system 
which is assumed uniform. 

Under hunting conditions Eq. 40 gives 


(42) 


— by p’ 


= At” + Adm: 


where 5,,, is defined in Appendix I. The torque equation gives 


Ar = yt" (At?) + dm At’ + (Sp? + Bo — (A). 


(43) 


In Appendix I, = 7% + + La, and = 72+ + La’, 
where 7” = 28 = i¥4/ v2, 1*4 denotes the direct-axis field excitation, and 
and 7° are the f, components. 


Part III 
6. THE QUASI-HOLONOMIC FRAME 


In the quasi-holonomic reference frame we shall keep the rotor 
currents and voltages and the mechanical coordinates (@ and w) as 
they are and transform the stator currents and voltages to a synchron- 
ously rotating reference frame. According to the revolving-field theory 
(3, 4, 5, 6), the transformed components are the forward and backward 
(f, 6) components (or Ku variables). According to the two-reaction 
theory developed by Blondel, Doherty, Nicle, and Park (10), the trans- 
formed components are the direct and quadrature (d, g) components 
(or Park variables). On the stationary reference frame, the phase 
voltages and currents of a three-phase synchronous machine can be 
transformed to their symmetrical (positive, negative, and zero sequence) 
components or to the alpha and beta (a, 8) components. (The in- 
stantaneous symmetrical (1,2) components are termed Lyon variables 
and the corresponding (a, 8) components are termed Clarke variables). 
For a two-phase machine, the armature phase currents and voltages can 
be labelled (a, 8) currents and voltages. The relation between the 
symmetrical (1, 2) components and the (a, 8) components is similar to 
the relation between the (f, 6) components and the (d, g) components. 

In the following sections we designate the indices a’’8’’ as the new 
indices for the quasi-holonomic reference frame similar to the use of 
a’8’ as the indices for the nonholonomic reference frame. Also we use 


the indices m’’n”’ to denote the electrical quantities and the indices u’’v 
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3 
: 
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j 
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: 
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to denote the mechanical quantities. We introduce the indices r’’s’’ to 
denote the rotor quantities alone. However, just as the indices a’’B”’ 
include both mn” and u’’v’’, it is understood that the indices m’’n”’ 
include the indices r’’s’’. In the quasi-holonomic reference frame here 
chosen, we let e,. =e,, = 7, = R,,, and Ly = L,,. 
Similarly, we let = = 2°, Jury = Ju, and By. = By. 

Let the transformation tensor be 


= f(x), = 6-7, and Cyr = bye (44) 


=z = and = = (45) 


By using double primes instead of single primes in Eq. 30 we get the 
equation of motion as 


_ at”, ar" 


Ox” 


+ Cy Car (= “oxy J (46) 


where 7”’ and F” are defined similar to Eqs. 31 and 32 with double 
primes replacing single primes. From an equation similar to Eq. 29 
we get the equation of voltage in a quasi-holonomic frame as 


din” 


In Appendix II it is shown that for the quasi-holonomic frame chosen 
= Ran, ANd Ann = Constant, and hence = 0. 
The equation of torque in a quasi-holonomic frame is given by 


(49) 


Note again that for the quasi-holonomic frame chosen R,»,., = Ru, 
= Aw, ANd Ann = Constant so that = 0. 


: 
; 

| hen 
: 
é 
: 

m 

at 
ox 
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: 
: 
ag 

kf’ 
n 0 x u’! 
: 
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In Appendix II we have shown that in terms of forward and back- 
ward (f, 6) components Eq. 47 gives 


di? diFé 
+ + Lyre dt ) 


dt 
+ + (50) 


= Ry + (1,2 di 


difé 


+ Lal + Lort®*)w (51) 


@ra = Rept + + Lr nt) (52) 
“a dt 


where w = 7“’’, and 
Ry; = Rwy = R,, = Lup = Ly, Lin Lu; 
Rrr R; Rra, Lrr = L; = Lra, 

Lyr = Lr; = = Lor = Lu. 


Note that in case we have more than one winding on the rotor, Eq. 47 
would contain the equation of voltage for the rotor as 


(53 
dt ) 


e 


In the case of a synchronous machine with a direct-axis field winding 
only, the impedance matrix is 


f b Fd 
f | (p + jw)L» | (p + jw)Ls 
| pi 


Ri + (p — jw)L2| (pb — jw)Ls| (54) 


R; + pL: 


= b jJw)L 


pls 


= d/dt, 


The equation of torque (49) becomes, for p 


(B + Jp)w + (jos) — 1°( joy) 

(B+ Ip)o + — i) + (itil — 
(B + Jp)w — — (Le — L,) ivi? 

= (B+ Jp)w — + 


+ 
ip 
| 
Fd | 
| 
a 
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where Loa = V2L4, (La — L,) = 2L2, = — — 2*)/v2, and 
+ 72°)/v2. Note that (La + L,) = where La and L, are 
the direct and quadrature axis inductances. The flux-linkages are 


defined as follows: 
Lake (56) 
= Ly? + Lol + (57) 
oa = Lat + 


@, = L*. 


7. A FOUR-DIMENSIONAL PROBLEM 


We shall consider Eqs. 50, 51, 52 and 55 as simultaneous equations 
of an electromechanical system with four coordinates x/, x°, x¥4¢, and @, 
and their velocities (or currents) 7/, 7°, 77“, and w. We have been able 
to reduce the angle-varying parameters to constant parameters. How- 
ever, if w is considered as a variable, Eqs. 50 and 51 contain terms that 
are products of 7/, 7°, or 2*¢, and w. Similarly Eq. 55 contains terms 
that are products of 7’, or 7°, and 7*“, and of (7’)? and (7°)?. We shall 
introduce the following notations : 
ghd and y= (60) 


v= 4 


fi = eb, fs = @€rad, and ts (61) 


The constants R,, R;, Li, Le, L3, and Ly are as defined before. We 
shall also keep the constants B and J for the mechanical coordinate 
equation. Thus we have 


dv! v? 
»= Le + Ly — + Rw? — j(Lw' + + Lw*)v' (63) 


dv! 


= dt 


La 


fs 


Since these differential equations are nonlinear, we cannot use the 
method of Laplace transforms or Heaviside’s operational calculus. 
Let p denote d/dt only. Then the above equations become, for w = pv, 


: (59) 

4 

: f 1 

: 

dv? dv3 
3 

fs L; — + Rw? (64) 

dt 

= J— + Bot t+ — + — (65) 

dt 
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+ Low? + + Ry! + j(Lw'! + Lev? + (66) 
= + + Lw*)v4 (67) 


fi 
fe 
= 
fi 


For a non-salient-pole synchronous machine, L. = 0, the above equa- 


= Lw'! + Lyw? + + Rw’? 


Lyw' + + + (68) 


+ Bot + — v*)v? + jL2(v'v! — (69) 


tions simplify to 


fi = Lyw'! + + Rw! + J(Lyw! + Lyv*)v' (70) 


fe = + + Ry! — + (71) 


fs = L + + + Rw’ (72) 


fs = + Bot + — 


Equations 66, 67, 68, and 69 can be transformed into the d, g com- 
ponents to give the following: 


= Lawt + Lraw?? + Ryv* — (74) 


L,wt + + (Laav?? + Lav*)v' (75) 


+ + R;vF4 (76) 
(La — L,)vw? (77) 


+ Bvt — 


where fa 44, pvt = wt, vt = 1%, put = wt, = = 


po¥4 = wF4, and Laa = V2Lu, Rea = Rz, and Lea = Ls. For a non- 
salient-pole machine, Za = L,, and the term for reluctance torque dis- 


appears in Eq. 77. 


8. THE TORQUE EQUATION 


We shall investigate the torque equation given by Eq. 55, or Eq. 
Let 7¥4 remain constant during the mechanical oscillation. Let 


77. 
also 


i4 = A cos8, —Asin# (78) 


where A is a constant and @ (or 6) is the angular displacement between 
a point on the shaft at any instant and a synchronously rotating 
This is known as the power angle. Equation 55 then gives 


frame. 


dt dt? 


“4 
‘ 
bs 
: 
3 
(73) 
A 
Ja 
3 
4 
: 
ae 
3 
is 
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where K = and = L.A? = A?(La — L,)/2. For a non- 
salient-pole machine, L. = 0, and hence = 0. Equation 79 is 
nonlinear on account of the presence of sin @ and sin 26 terms. For 
K, = 0, Eq. 79 is similar to the equation for a pendulum with damping 
constant B and spring (stiffness) constant K. A more elaborate form 
of Eq. 79 was given as (11, 12, 13, 14, 15) 


120 
p= + B(1 — bcos 26) + 20. (80) 


As Eq. 80 and a simplified form of Eq. 79, with A, = 0, were treated in 
(9), we shall discuss the solution of Eq. 79 here. The constants J 
and K can be normalized to unity, after proper changes in the other 
constants and the time scale and the torque unit. The new equation 
becomes 


+ + sin + Ky’ sin TFT. (81) 
1 


We have retained the symbol 7 for torque in spite of the change of unit. 
Note that 7» denotes the original torque load, and 7, denotes the addi- 
tional torque suddenly applied. The phase-plane equation for the 
phase plane (8, w) is 
dw 7 — sin 26 — sin@ — B’w 
(82) 
w 


To find the saddle point, set w’( = dw/dt) and w to zero. Thus 
7 = K,’ sin 26, + sin 6, (83) 


where 6, denotes the angle at which the saddle point is located. The 
slope at the saddle point is given by 


S = —a+t — 2K,’ cos 26, — cos (84) 


where a = B’/2. If Eq. 80 is used as the starting point, 
a = B’(1 — bcos 26,)/2. If @ is expressed in degrees, the terms 
(— 2K,’ cos 20, — cos @,) must be divided by 57.3. For illustration, 
let B’ = 0.01, K.’ = 0.35, and + = 1.0. Then Eq. 83 becomes 
1 = 0.35 sin 26, + sin 6, and 6, is obtained as 90°. The slope at the 
saddle point is 


S = — 0.005 + [(0.005)? + 0.70/57.3 2 
= — 0,005 + 0.110 = 0.105 or — 0.115. (85) 


| 
| 
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Following the negative slope, the w — 6 trajectory is traced backward. 
This gives the separatrix which separates the stable region from the 
unstable region. The area inside this separatrix is stable. Beyond the 
separatrix the system is unstable and the machine pulls out. The 
construction of the trajectory can be made either by the zsocline method 
of the acceleration-plane method (11, 12). Suppose that the separatrix 
trajectory reaches w = 0 at 6 = 4. Then the initial load is* 


To = 0.35 sin 26, +. sin Ao. (S6) 


The limiting value of the additional torque load is 


= 1 — To. (87) 


9. THE VOLTAGE EQUATION 


Now if we assume that w is constant in Eqs. 50, 51, and 52, the 
impedance matrix given in Eq. 54 can be utilized to write the following 
equations : 


[Ri + (p + jo)Li + (p + + (p + (88) 
(p — + [Ri + (p — jw)Li + (p — (89) 
(90) 


Since the above three equations are linear, superposition holds. We 
may put éra = 0 and study the effect of suddenly applying e; and e, on 
7’, 7*, and the transient component of 7*¢. From Eq. 90, for era = 0, 
we get, after transposing terms, 


ey 


Crd 


+ pLs)t?4 = — pL’ — pla’. (91) 


Premultiplying Eq. 91 by pL, gives 


pL.(R; + = — + 2). (92) 


Since w is considered constant, premultiplying Eqs. 88 and 89 by 
(R; + pL;) and substituting into the new equations Eqs. 91 and 92 give 


(R3+pL;)e;= (Ri +pLs)[Rit (p+ jw)Li]—p(pt+ jw) Le jv 
+ (93) 


(Ri; +pLs3)e,= (p— jo)Li]—p(p— jwo)Le 
jw) (94) 


31f ro is given as 0.434 per unit, Eq. 86 gives 6) as 15°. For a total torque load of say 
r = 0.95 corresponding to 7; = 0.516, the phase-plane trajectory can be constructed according 
to Eq. 82. Then w(t) and @(t) can be obtained from the trajectory by the relation (A@/At) = war, 
where w,, denotes the average angular velocity during the small interval. 


= 
= 
|_| 
| 
; 


We have eliminated one column and one row from the matrix (54) in- 
volving the field quantities and obtained the operational (or transform) 
impedance matrix 

b 


f | Rit (p+ jo)Lil(p) | (pb + jw)L2(p) 


b | (p — jw)L.(p) R, + (p — jw)Li(p) 


where 
2 
pe se (96) 


Li (p) = R; rn pL;’ L2(p) R; +4 pL; 


For getting the transient inductance (or reactance), p is set to ~. 


Hence 
Lé ) Q” 


The sum of L,(p) and L2(p) gives La(p), and the difference between 
Li(p) and Le(p) gives L,(p). Thus 


For p = ~%, we get 


La 
= kL L,’ = L L,’ L, 99 
1 + Le LaL ra 


where Lrg = Ls and Lua = V2L,. If there are damper windings, L;(p) 
and L2(p) will have more complicated expressions, and so will be La(p) 
and L,(p). The subtransient inductances are again obtained by setting 
pto x. Thesymbols used are L.’’, La’, and L,”’ (see Appendix I). 

In terms of d, g components with field coordinate eliminated, we have 


d q 


d| Ri, + pla(p) — wL,(p) 


wLa(p) R, + pL, (p) 


The transient analysis of synchronous machines can then be made with 
the help of either Eq. 95 or Eq. 100. 


: 
ate 
(95) 
| 

pL ad” 
: 
La(p) = La L,(p) = Ly. (98) 
R 3 pL 3 
: 

: 
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CONCLUSIONS 


An electromechanical system can be expressed by Lagrange’s equa- 
tion of motion in the holonomic frame. As the angle @ is the ‘‘fourth 
dimension,” inductances which vary with @ have to be considered 
nonlinear. By transforming the system from the holonomic frame to 
the nonholonomic frame, the new equation of motion takes a more 
general form involving the presence of the affine connection which re- 
duces to the Christoffel symbol in the holonomic frame. ‘Theory of 
hunting can be developed based on the nonholonomic frame. If we 
keep the rotor quantities and the mechanical coordinate as they are, 
we get the quast-holonomic frame. Then certain nonholonomic objects 
become tensors. By referring the stator quantities to the rotor, we 
have two sets of components: the f, ) components (or Au variables) and 
the d, g components (or Park variables). The equations of voltage 
involve the products of 7’, 7, or 2*4, and the angular velocity w. Since 
we consider w as a dependent variable in the “fourth equation” (see 
qs. 66-69), these equations 66 and 67 are nonlinear. The equation 
of torque involves the product of 7¢ and 7¥¢ in the case of a non-salient- 
pole machine and of 74 and 7¢ additionally in the case of a salient-pole 
machine. Thus Eq. 69 is nonlinear. (Only Eq. 68 is linear in form). 

The equation of motion of an electromechanical system is in general 
nonlinear either in the holonomic frame, or in the quasi-holonomic frame, 
or in the nonholonomic frame. Consider the set of four equations 
(66-69) again for the quasi-holonomic frame, we find that the torque 
equation (69) alone gives the well-known pendulum equation (79) for a 
non-salient-pole machine. The solution previously obtained from a 
differential analyzer can now be obtained by the phase-plane method 
used in nonlinear analysis. The effect of saliency can be taken into 
account by the form shown in Eq. 80. Stability study can be made by 
finding the saddle point and constructing the separatrix. If one is only 
interested in electrical transients under the assumption that w is con- 
stant, then Eqs. 66-68 form a set of three /inear differential equations 
with vt = w considered as a constant. Applying the principle of super- 
position and eliminating 7¥¢ = v’, we get an operational (or transform) 
impedance matrix given by either Eq. 95 or Eq. 100. Transient an- 
alysis can thus be made by transforming the equations of voltage from 
the holonomic frame to the quasi-holonomic frame with great advantage. 
It is suggested that Eqs. 66-69 or 74-77 should be solved simultaneously 
for the general case where separate solutions obtained by dividing the 
four dimensions into two sets (3 and 1) are no longer valid or justified. 
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APPENDIX I 


Starting from the stationary frame, the armature (stator) phase equations can be given 


as follows: 


Ry + poi (AI-1) 
Ri? + poe Cs (AI-2) 


where 7! and @ denote the positive and negative sequence armature current components, ¢; 
and ¢2 denote the positive and negative sequence flux-linkages, and e; and ez denote the positive 
and negative sequence voltage (terminal) components. The flux-linkages are (5): 


= + Lit + Lote? (AI-3) 
= + Lit? + Lette (AI-4) 


where 6 = wt, L; denotes the average value of La and Lg, the direct and quadrature inductances 
of a salient-pole machine, L2 denotes one half of the difference (La — Lg), and Lea is the mutual 
inductance between the armature phase and the direct-axis field. The field current com- 
ponents are defined by 


(AI-5) 
where 1*¢ and 7"? denote the field excitation currents in the direct and quadrature axes. The 
field flux-linkage components are similarly 


v2 F v2 


Since 1¥@ is zero, Eqs. Al-5 and AI-6 simplify to 


= 18 = or = = Laat* = 


In the literature, the value of Lag is often assigned as unity in the per-unit system. 
Transforming Eqs. AI-1 and AI-2 into a rotating reference frame gives 


Ri! + poy + =e (AI-8) 
Rit? + pdr — (AI-9) 


oy = or + Li Lat? (AI-10) 
on = + Lot! + (AT-11) 
(AI-12) 


(Al-6) 
where 
4 
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er = = (AI-13) 


(AI-14) 


os = dy = gre". 


Note that for irg = 0, and Laa = 1, dr and ¢p are defined by Eqs. (AI-7) and may be replaced 


by 7” and 74. 
Assuming all field axes eliminated, the conventional equations of hunting are given as* 


+ (p + jw) + + —-(AT-15) 
+ (p — jw) + — (0). (AI-16) 


ber — (bp + jw) Adr 


ben — — jw) Ads 


From the torque equation 


= (Jp? + Bp)o — i? jby + = (Jp? + Bp)O — i*By — By 


(AI-17) 


where J denotes moment of inertia, B denotes friction constant, By = j@y and By = — jdp, 


we get the equation of hunting torque 


(Jp? + Bp)(40) — jt’(Adr) + ji! (Ads) + + (AI-18) 


Ar 


where 


= — — + (AI-19) 


by = — By + — 


by = — By + — = — + (AI-20) 


Transposing terms in Eq. AI-18 gives 


Ar + ji>(Adr) — = + + (Jp? + Bp) (d8). 


(AI-21) 


Note that the second and third terms on the left side of Eq. AI-21 are equivalent to 7%Atra, 
where 14 = (— ji/ + jz®)/v2, and Agr = Ads = At¥4/v2 for Laa = 1. 
The absolute differentials are related to the conventional changes as follows: 


= Ail — jit (A@) (AI-22) 


51> = Ar? + (AI-23) 


Similarly we have 


Bey = Aey — (AI-24) 


de, = + Ad). (AI-25) 


Substituting Eqs. AI-22—AI-25 into Eqs. AI-15 and AI-16 gives 


Aes — (p + jw) At¥4/y2 
= [Ri + (pb + jo)Li" Jai! + (p + jw) At? — (p + (AI-26) 


— (p — jw) Ai¥4/v2 
= [Ri + (p — Jai? + (p — jw) Ail — (p — (AL-27) 


where by and }, are defined in Eqs. AI-19 and AI-20. Similarly, Eq. AI-21 gives 


Ar + = + + (Jp? + Bp — jbvif + jbyi?) (a8). (AI-28) 


Note that in using Eqs. AI-24 and AI-25, es should be replaced by (Rif + jwos) = (Rit! + wBy) 
and e, should be replaced by (Rit? — jwds) = (Rit? + wBy). The incremental field excitation 
current A7¥¢ represents Aera/Rra, where era denotes the direct-axis field excitation voltage and 
Rra denotes the field winding resistance. The effect of the field winding and the damper wind- 
ings is taken account of in the use of “transient’”’ or “‘subtransient’’ inductance ZL," and L2” 


corresponding to La” and L,”. 


* $i and ée are the absolute differentials and A@ denotes the conventional change of angle. 
The symbols LZ,” and L2"’ are introduced to denote the ‘‘subtransient”’ inductances correspond- 


ing to L,’’ and L,”’ used in two-reaction theory. 
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APPENDIX II 


Equations of Voltage for Transient Analysis 


From Eqs. AI-1—Al-4 we get 
= Ry + p( Lit! + Leet? + Lyeiti?4) (AIT-1) 
= Ry + p(Lit® + Loe + Lye 4) (AIT-2) 
where Ls = Lea/V2, and i*4 = y2i* = y2i%. If we assign the index 3 to the field excitation 


current 7*4 in addition to the indices 1 and 2 for the positive and negative sequence com- 
ponents, we have 


Ri, = Ro = R,, Li = Lo» = Ly (AIT-3) 
Lie = Lay = Loe 20 Ly = Lye”, = Liye 19 (ATT-4) 
Note that Ly», is a function of @for m # n, and is nota function of the current 7". ‘The direct- 


axis field or rotor equation is given by 
= era = Ry + + Lee! + Lee? (ATT-5) 
where 7? = i¥4 and e; = era. Thus we have 
R33 = R; = = = Ly (AII-6) 


Ls 


L.e~**, = Lye**, (AIT-7) 
Note that Z,; and L3; are conjugates, and Le; and L32 are conjugates. Similarly, Lis and Lo 
are conjugates. 

Let the symmetrical components be transformed to the f, 6 components by 


Ci= |“ (AII-8) 


lhe field current component will remain the same. We shall use the index 3 for the original 
frame and the index F for the new frame in combination with the f, 6 components. 


= 3 1 (AIT-9) 
Then we have 
F 
1 | 
3 


Let the positive and negative sequence flux-linkages be defined by 
= Ly + + Lyeiti? = (AIT-11) 


o2 = Ly? + Loe + Lye it = (AIT-12) 
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Then it is seen that by Eqs. AII-8 and AIT-9 or by Eq. AIT-10 we get 


oy = Lyil + Loih + LyiF4 (AIT-13) 


+ Lol + LaF 4. (AIT-14) 


Thus we have 


(AIT-15) 


and Lyr 


+ Loy = Le, 


Ly = Lw = Li, 


Transforming Eqs. AII-1 and AII-2 to the rotating reference frame (5) gives® 


ey = Ry! + (p + jwldy = Ri! + poy + wBy (AII-16) 


ep = Ry’ + (p — = Rit” + + (AIT 17) 
where 
By = joy, By, = — joo. (AIT-18) 
Note that 
Ry Ry @ (AII-19) 


The field equation (AII-5) becomes 


era = + + Lal + Lai’). 


(AIT-20) 


Thus we have 


Rrr = R;, Lrr = Ls, = Ley = (AI[-21) 


So by choosing a quasi-holonomic frame we have succeeded in transforming the equations of 
voltage to a new set of equations with constant as well as reciprocal coefficients. Whether 
these equations are linear or nonlinear, however, depends upon whether the angular velocity 
w is considered as another dependent variable (in the fourth dimension) or a mere constant. 

Equations AII-16 and AII-17 can be rewritten as follows: 


(AII-22) 
(AIT-23) 


Ry (p + jw) (Lyi! + Lot + 


= Ry? + (p — jw) (Lit? + Lai! + 


In terms of d, g components these become 


(AII-24) 
(AIT-25) 


Riit + poa + wBa 


ea = Rit + p(Lait + Laat? 4) — 


ey = + p(Lgit) + w(Lait + 4) 


where 


oa = Lait + Leai®4, = Legit (AII-26) 


and 


Ba B, da. (AIT-27) 


The field equation can then be rewritten as 


eva = + + Laat). 


(AIT-28) 


If the angular velocity w is assumed constant, the effect of any disturbance in the terminal 
voltages calling for a transient analysis of the electrical quantities can be studied by the equa- 
tions of voltage given in terms of the f, 6 or d, g components on a rotating reference frame. 
The insertion of external networks and transmission lines will modify these equations. How- 
ever, if these external networks and transmission lines are also referred to the rotating reference 
frame, and represented by A and B, then Eqs. AII-22 and AII-23 can be modified by adding 


A and B to R,, respectively. 


5 It is easy to verify that 
= pos + = + jw)dy 
p( gre = — joww) = — 


where w = 


~ 
= Lip = — 
| 
ber 

: 

Ss 

i 

2 


H. Ku 


REFERENCES 


(1) G. Kron, ‘“Quasi-Holonomic Dynamical Systems,"’ Physics, Vol. 7, pp. 143-152 (1936) ; 
“A Short Course in Tensor Analysis for Electrical Engineers,’"’ New York, John Wiley 
& Sons, Inc., 1942; ‘‘Tensors for Circuits,"’ New York, Dover Publications, Inc., 1959. 

(2) L. V. BEw ey, “Tensor Analysis of Electric Circuits and Machines,’’ New York, The 
Ronald Press Company, 1961. 

(3) Y. H. Ku, “Transient Analysis of A-C Machinery,” Trans. ATEE, Vol. 48, pp. 707-715 
(1929). 

(4) Y. H. Ku, “Transient Analysis of Rotating Machines and Stationary Networks by 
Means of Rotating Reference Frames,’ Trans. ATEE, Vol. 70, Part 1, pp. 943-957 
(1951). 

(5) Y. H. Ku, “Electric Energy Conversion,’’ New York, The Ronald Press Company, 1959. 

(6) Y. H. Ku, “‘Rotating-Field Theory and General Analysis of Synchronous and Induction 
Machines,”’ Proc. JEE (London), Vol. 99, Part IV, pp. 410-428 (1952). 

(7) G. Kron, “Equivalent Circuits for the Hunting of Electrical Machinery,"’ Trans. AITEE, 
Vol. 61, pp. 290-296 (1942). 

(8) G. Kron, “A New Theory of Hunting,’”’ Trans. AJTEE, Vol. 71, Part III, pp. 859-862 
(1952). 

(9) D. C. H. H. Woopnson, “‘Electromechanical Energy Conversion,’"’ New York, 
John Wiley & Sons, Inc., 1959. 

(10) R. H. Park, “‘Two-Reaction Theory of Synchronous Machines—I and II,” Trans. AIEE, 
Vol. 48, pp. 716-729 (1929); Vol. 52, pp. 352-355 (1933) 

(11) Y. H. Ku, “Nonlinear Analysis of Electromechanical Problems,”’ Jour. FRANKLIN INST., 
Vol. 255, pp. 9-31 (1953). 

(12) Y. H. Ku, ‘Analysis and Control of Nonlinear Systems,’’ New York, The Ronald Press 
Company, 1958, pp. 108-112. 

(13) W. V. Lyon anp H. E. EpGerton, ‘Transient Torque-Angle Characteristics of Syn- 
chronous Motors," Trans. ATEE, Vol. 49, pp. 686-698 (1930). 

(14) H. E. EpGerton anp P. FourMarter, “Pulling Into Step of Salient-Pole Synchronous 
Motor,” Trans. AIEE, Vol. 50, pp. 769-781 (1931). 

(15) D. R. SHoucts, S. B. Crary anp A. H. Lauper, “Pull-In Characteristics of Synchronous 
Motors,” Trans. AJEE, Vol. 54, pp. 1385-1395 (1935). 

(16) Y. H. Ku, “Transient Circuit Analysis,’ Princeton, D. Van Nostrand Company, Inc., 

1961. 


: 
4 
: 
: 
: 4 


infinite metal target. 
connection with the shaped-charge effect (1),? the impact of cylinders 


PENETRATION OF A ROD INTO A SEMI-INFINITE TARGET 
BY 


WILLIAM A. ALLEN! AND JAMES W. ROGERS! 


ABSTRACT 


The penetration of metal rods into semi-infinite metal targets has been investi- 
gated experimentally at velocities up to 0.3 cm/ysec. The rods were composed of 
Au, Pb, Cu, Sn, Al, and Mg; the targets were aluminum. Results are compared with 
predictions from the hydrodynamic theory of jet penetration. Basic assumptions of 
the hydrodynamic theory were used to determine an effective yield strength of the 
target-rod combination. Experimental results indicate that the effective yield 
strength is relatively independent of the strength of the rod. The hydrodynamic 
theory of penetration was determined to be generally acceptable except where the 
density of the jet is much greater than that of the target. A gold jet reveals a new 
effect of secondary penetration which results in penetration greater than that pre- 


dicted by theory. 


INTRODUCTION 


This paper deals with the penetration of a metal rod into a semi- 
The problem has been considered previously in 


against plates (2), and the mechanism of indentation in a hardness 


test (3). 


Impact velocities produced by steel shaped charges ranged 


up to 0.8 cm/ysec, the cylinders (2) were fired against plates at veloci- 
ties up to 0.02 cm/yusec, and the hardness tests were conducted at 
essentially static conditions. A hydrodynamic theory has been intro- 
duced to explain the penetrations at shaped-charge velocities; the 
theories of elasticity and plasticity have been used in the low velocity 


range. 
the validity of the hydrodynamic theory of penetration. 


The experiment discussed in this paper was designed to check 


HYDRODYNAMIC THEORY OF PENETRATION 


length /, velocity v, and density p, into a semi-infinite target of density p» 


can be readily deduced. 


density p, traveling with velocity v toward a semi-infinite target of 


density 


Figure 1c indicates the actual penetration process where the bottom of 


the crater travels at velocity wu. 


of primary penetration; the jet has been converted entirely into a 
cylindrical shell which, in general, is traveling with respect to the target. 


If the movement of the jet material is toward the bottom of the crater, 


1U.S. Naval Ordnance Test Station, China Lake, Calif. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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The conventional theory (1, 4, 5, 6) of the penetration of a jet of 


Figure la represents a rod of length / and 
Figure 10 illustrates the situation just after impact. 


Figure 1d represents the termination 


: 
= 
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the possibility exists that subsequent secondary penetration will occur. 
Assume that penetration begins at the initial impact, that the penetra- 
tion rate u is constant, and that penetration ceases when the rod is 
used up. If initial transients, secondary penetration, and cavitation 
are neglected, the duration of the penetration event is given by 


l 
(1) 


t = 
(v — un)’ 


the final crater depth P = wt is given by 


(d) 


Fic. 1. Hydrodynamic penetration of a rod into a semi-infinite target. 


and the stagnation stress is given by the modified Bernoulli’s equation 
spilv — = dpou? + ¢, (3) 


where the term ¢ can be considered a dynamic yield strength of a solid 
target with respect to a fluid jet. The quantity @ is regarded as a 
strain-rate function involving the material strength of both rod and 
target. However, if ¢, determined experimentally, should turn out to 
be relatively independent of jet material, then it can be inferred that 
the strength of the rod may be neglected and the resultant value of ¢ 
is a material characteristic of the target material alone. Equation 3 
is valid for incompressible steady-state flow, but can be applied to 
cases of compressible isentropic steady-state flow. For the case of 
subsonic flow of a jet, Eq. 1 is a good approximation but is not exact 


2 

: 

ul 

(v — u) 

Uf 

Y Yy 

= P. Yi 
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because the particle velocity component parallel to the trajectory is v 
only at infinity. 
Solve Eq. 3 to obtain 


pw — — 2(p1 — — @) }'? 
Pi — Pe 


“u = 


Substitute Eq. 4 into Eq. 2 to obtain 


P=! 


In the limit of high velocity the effect of ¢ becomes negligible and Eqs. 4 
and 5 approach the respective forms 


“= (6 
1 + (pi/p2)'” 


P = (7) 


The quantity ¢@ in Eq. 3 can be defined for the critical case where 
u = 0 by means of the relation @9 = $piv)?._ The critical velocity v» is 
the striking velocity that must be achieved before steady-state pene- 
tration can occur. 

The phenomenon of jet penetration is related to the mechanism of 
indentation in a hardness test (3, 7). A paper has appeared (8) in 
which an approximate solution was given for a rigid cylindrical punch 
penetrating an ideal plastic solid. Plastic yielding occurred when the 
mean pressure on the punch is about 2.8 Y where Y is the yield stress of 
the material. A more recent analysis (9) for the case of a spherical 
punch penetrating an ideal plastic solid indicates that @» = 2.6Y. 
The static yield strength of 7075—T6 aluminum (10) is 0.00496 megabar. 
Data have been published (11) which indicate that the ultimate dynamic 
yield strength of duralumin of Brinell hardness 133 is 50 per cent higher 
than its static value. Thus, a theoretical value of @) can be obtained. 
For future reference assume ¢o = 0.019 megabar for 7075—T6 aluminum. 

The conventional hydrodynamic theory of a rod penetrating a 
target has never been adequately verified by shaped-charge jets because 
the assumed ideal conditions of the theory are difficult to achieve experi- 
mentally. A shaped-charge jet does not have a uniform cross section. 
In fact, the jet often consists of a stream of discrete fragments. Even 
when a jet remains continuous a velocity gradient exists along its 
length. Both the length / and velocity v of a shaped-charge jet are 
rariables. Consequently, accurate experimental verification of Eq. 5 
has not been reported previously. 


| 
| 
vil — = — 1/2 
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EXPERIMENTAL RESULTS 


A paper has appeared (12) that describes development of a hyper- 
velocity precision impact instrument. This research tool has been 
used to check the theoretical penetration relation expressed by Eq. 5. 
A short metal wire segment approximately 0.020 in. in diameter was 
aligned on the trajectory of the instrument by means of a precision jig. 
This wire segment will subsequently be referred to as the jet. A 
7075—T6 flat-nosed aluminum projectile, subsequently to be referred 
to as the target, was fired at the end of the jet (13). After impact the 
target travels down the range and is captured without appreciable 
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Fic. 2. Penetration of jets of different metals into 7075-T6 aluminum. 
The curves are empirical. 


deformation in a cotton backstop. The target is then sectioned in order 
to determine the crater volume and the penetration. A cross section 
of a typical crater produced by a lead jet appears as Fig. 6 in (12). 

Figure 2 displays experimental data for penetration of jets of various 
metals into 7075-T6 aluminum. The ordinate of Fig. 2 is relative 
penetration y defined by the relation 


(8) 


where P is the measured penetration, / is the length of the jet, and 
Pi, p2 are the respective densities of jet and target, The solid curves 


: 
C25 
: 
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in Fig. 2 are empirical relations of the form 


(2) (1 a) (v Vo) + Vo (9) 


where a and v» are adjustable parameters. <A plot of Eq. 9 for the case 
of the lead jet appears as Fig. 7 in (12). The intercept vo on the x axis 
for each curve is a threshold velocity that must be achieved before 
steady-state penetration can occur. To each value of vo is associated 
the critical stagnation stress @o = 3pit0’. 

Table I includes the physical properties of the jets used in the experi- 


TABLE |.—Physical and Experimental Properties of Jets. 


Critical 

Length? /, Diam d, Density 1, velocity v9, Adjustable 
Metal cm cm g/cc cm/psec parameter a 
Au 0.467 0.0635 19.30 0.0479 0.90052 
Pb 0.465 0.0737 11.09 0.0583 0.81329 
Cu 0.432 0.0635 8.96 0.0672 0.81705 
Sn 0.540 0.0711 7.30 0.0675 0.76894 
Al 0.801 0.0787 2.70 0.1170 0.68072 
Mg 0.796 0.0635 1.74 0.1481 0.64308 
@ Nominal value. The length of the jet varied slightly for different rounds. : 


ments. The last two columns of Table I present the adjustable param- 
eters Yo and a that appear in the empirical curves of Eq. 9. The ex- 
perimental data of Fig. 2 can be displayed in an alternative and sugges- 
tive manner. Equation 2 written in the form 


Pv 


) 
(P +1) 


ut 


can be used to calculate the penetration velocity « as a function of 
the striking velocity v. Eliminate y between Eqs. 8 and 9 and substi- 
tute into Eq. 10 to obtain the relation 


u = a(v — (11) 


Figure 3 is a plot of Eq. 11 for the four jet materials of Fig. 2. The 
values of the adjustable parameters a and vp in Table I were determined 
originally by fitting a straight line by the method of least squares to 
data plotted in the form of Fig. 3. 

Substitute Eq. 11 into Eq. 3 and solve to obtain the relation 


1 1} 
a a 
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CM/p SEC 
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PENETRATION VELOCITY u, 
ro) 


| 
0.1 0.2 
PROJECTILE VELOCITY v, CM/uSEC 
Fic. 3. Penetration velocity as function of projectile velocity for jets 
into 7075-T6 aluminum. 


Too much physical significance should not be ascribed to Eq. 12 since 
its functional form is a simple consequence of the linear relation ex- 
pressed by Eq. 11 and plotted in Fig. 3. Write Eq. 12 in the form 


o(u) = do + bu + cx, (13) 


where the constants ¢o, 0, and c were calculated from Eq. 12 and are 
listed in Table II. Figure 4 is a plot of Eq. 13 where the values @o, 0}, 


TABLE II.—Parameters that Specity Dynamic Yield Strength in Eq. 13. 
Jet Material do b c 
Gold 0.02212 0.1022 —1.27¢ 

Lead 0.0189 0.1487 —1.10 
Copper 0.0202 0.1348 —1.17 

Tin 0.0166 0.1481 — 1.06 
Aluminum 0.0185 0.1482 —1.09 
Magnesium 0.0191 0.1430 —1.12 
Average 0.0187 +0.0012 0.144+0.005 —1.11+0.03 


® Rejection"of these data in averages was justified by Chauvenet's criterion. 


c and the indicated dispersions were obtained from Table II. The 
average value ¢) = 0.0187 + 0.0012 megabar in Table II can be com- 
pared with the theoretical value ¢) = 0.019 megabar for 7075—T6 
aluminum. The quantity = 0.144 + 0.005 can be considered a 
measure of the strain-rate effect. The quantity c = — 1.11 + 0.03 
substituted into Eq. 13 forces the curve to reach a maximum value 


2 
280 
OAu DAI 
9 O O 
e 
Mg 2 
: O O O 
O O O 
O 
: 
: 
; 
- 
: 
= 


Oct., 1961.] PENETRATION INTO A SEMI-INFINITE TARGET 


- mM 


Static yield stress 


— 


DYNAMIC YIELD STRENGTH ¢, KILOBARS 


fe) | | | | 
0.00 0.02 0.04 0.06 0.08 0.10 
PENETRATION VELOCITY u, CM/xzSEC 


Dynamic yield strength ¢ as function of penetration velocity for 7075-T6 aluminum 
target. Six jet materials were used. 


@ = 1.2¢0 at about u = 0.06 cm/ysec. The curve in Fig. 4 was dis- 
continued at the penetration velocity « = 0.10 cm/ysec because data 
for aluminum and magnesium have not been obtained for velocities 
in excess of this limit and the data for gold beyond this value were 
suspect. 


COMPARISON OF PENETRATION RESULTS WITH HYDRODYNAMIC THEORY 


Agreement between the experimental data and hydrodynamic theory 
will now be investigated. In view of the relatively large dispersion in the 
curve of Fig. 4, a tentative average value of ¢ for all metals is assumed 
based on the relation 


(14) 


Evaluation of Eq. 14 using the information of Table II leads to the 
value (¢) = 0.0222 megabar. Substitute this value in Eqs. 4 and 5 
and evaluate to obtain the curves of Figs. 5 and 6. The average value 
@ = (¢) implies existence of an average threshold velocity for each 
metal specified by the relation 


2P1(Vo)”. (15) 
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Both Fig. 5 and Fig. 6 reveal systematic deviations from hydrodynamic 
theory in the case of a high density jet such as gold. The actual pene- 
tration velocity « appears to be more accurately represented by the 
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Fic. 5. Comparison of experimental penetration data with hydrodynamic theory. The curves 
are specified by Eq. 5 where ¢ = (@) is the same constant for all metals. 
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Fic. 6. Comparison of experimental penetration data with penetration theory. The curves 
are specified by Eq. 4 where @ = (¢) is the same constant for all metals. 
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straight lines of Fig. 3 than by the curves of Fig. 6. It should be re- 
called, however, that the four curves of Fig. 3 were specified by eight 
adjustable parameters while the curves of Fig. 6 were fixed by a single 
constant (¢) which can be determined theoretically. 


ANOMALOUS PENETRATION OF A GOLD JET 


The anomalous penetration achieved in 7075—T6 aluminum with a 
high density jet such as gold can be explained by a new phenomenon 
of secondary or residual penetration. ‘This effect is essentially different 
from a phenomenon previously reported (14) under the same name. 
This previous effect was explained by target flow subsequent to jet 
penetration ; the discovered effect is a consequence of the difference in 
densities between the jet and the target. 

With respect to a coordinate system traveling with the velocity u 
of the bottom of the crater the jet velocity is v — u. The speed of a 
free streamline is not affected by a change in direction, and thus the 
velocity of the jet reflected from the bottom of the crater is given by 
u —v. In the coordinate system of the target the velocity of the 
reflected jet is 2“ — v. When the jet has used up its entire length / 
its mass retains this velocity. The densities of jet and target are usually 
comparable in shaped-charge experiments and the residual velocity 
2u — v is negligible compred with v. In the case of high density jets, 
however, the residual velocity is appreciable. The residual velocity 
v, = 2u —v can be calculated from Eq. 4, and is specified by the 
relation 


_ 2[ — 2(p2 — pi) (@ — — (02 + (16) 
— 


The threshold velocity vo is given by Eq. 15 where the brackets have 
been omitted without ambiguity. Ifv, > vo, secondary penetration can 
occur with magnitude given by Eq. 4 where v, calculated from Eq. 16 
is substituted for v. Equate Eq. 16 to v, eliminate ¢ = (¢) by means 
of Eq. 15, and solve to obtain 


Vo(3p1 + (17) 
pi 


Equation 17 specifies the critical velocity v. at which secondary penetra- 
tion begins to occur. This velocity occurs for the case of a gold jet 
into 7075-T6 aluminum at the value v, = 0.175 cm/sec; the curve in 
Fig. 7 illustrates the effect of secondary penetration in this case. The 
gold curves of Figs. 3 and 6 have not been extended past the critical 
velocity v,. A gold jet is the only metal that causes secondary penetra- 
tion in the velocity range studied, 
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Fic. 7. Penetration of a gold jet into 7075-T6 aluminum. The line is a theoretical curve that 
illustrates the effect of secondary penetration. 


Evidence exists that the anomalous penetration is not a result of 
cavitation. Penetration that results from cavitation would be inde- 
pendent of the length of the jet ; its value is sometimes taken equal to the 
radius of the cylindrical crater (14). The new effect is proportional to 
jet length. 
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TWO ANALOGIES IN SOLID CONTINUUM MECHANICS 
BY 
S. F. BORG! 


SUMMARY? 


This paper shows that the equations governing the behavior of moment and torque 
loaded thin plates which have small deflections are analogous to the equations of the 
two-dimensional linearized mathematical theory of isotropic elasticity. Because of 
the analogy, all solutions of two-dimensional plane-strain or plane-stress problems 
may be taken as solutions of a particular thin plate problem and vice versa. 

A three-dimensional generalization of the above is also developed. In this form 
the theory draws an analogy between solutions in three-dimensional elasticity and 
those of thin plates subjected to lateral and shear loads as well as moments and 
torques. 

The latter development is based on a tensoral generalization of the two-dimen- 
sional result to three dimensions. Beczuse of this analogy all the known solutions 
for thin plates subjected to lateral loads represent possible solutions in a three- 
dimensional elastic structure,’ and vice versa. 

Examples of both analogies are given. 


THE EQUATIONS OF THE TWO-DIMENSIONAL LINEARIZED MATHEMATICAL 
THEORY OF ISOTROPIC ELASTICITY 


The fundamental tensors of two-dimensional linearized elasticity 
are the stress tensor, T, given by‘ 


r-(% 
Tyx Gy 
and the strain tensor, 9 given by 
au 2+ 2)| 
oy Ox 


2\ ax ay F) y 


Note that these tensors are symmetrical. Also they are tensors in 
the sense that they satisfy a particular rotation of axes law, which 


1 Head, Civil Engineering Department, Stevens Institute of Technology, Hoboken, N. J. 

2 A preliminary form of the theory presented in this paper was given some years ago in 
the author's book ‘An Introduction to Matrix-Tensor Methods in Theoretical and Applied 
Mechanics,” Ann Arbor, J. W. Edwards, 1956. 

3 An analogy between the plate solutions and membrane solutions is discussed at some 
length in S. TrmosHENKO, “Theory of Plates and Shells,’’ New York, McGraw-Hill Book Co., 
Inc., 1940. 

4 A list of the symbols used in this paper appears in the Appendix. 
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typically, is 


T’ = RTR* (3) 


in which 


T is the tensor referred to the x, y axes 
T’ is the tensor referred to the x’, y’ axes (that is, the axes rotated 
about the origin) 

R is the rotation matrix defined by 


4) 


R* is the adjoint matrix, that is, the matrix R with rows and columns 
interchanged, and 
l;; is the cosine of the angle between the 7’ and the 7 axes, where 
i’ = x’, y’ and j = x,y 
The fundamental equations governing the behavior of a deformable 
body, in the two-dimensional linearized theory of elasticity, are: 


stress equilibrium, body forces neglected 


div T = 


z z 
Ox Tex: We 


boundary conditions on stresses 


or 


TN (7) 


og T l x 

tin m y 
in which (/, m) are the direction cosines of a normal to the boundary 
and X, y, are the boundary stresses in the x, y directions. 


or 


Hooke’s Law in its linearized isotropic form 


anda 


strain compatibility equation, in its linearized form 


(J. F. 1. 
: 
: 
: 
; 
(5) 
E E 2 
2 2 9 
oy? 
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In Eq. 9, v is Poisson’s ratio, E is the isotropic modulus of elasticity, 
J, is the trace of the stress tensor and & is the two-by-two unit matrix. 


THE FUNDAMENTAL EQUATIONS OF SMALL DEFORMATION THIN PLATE THEORY 


We may represent the loading, moments and torques acting on a 
typical element of a thin plate asin Fig. 1. The moments, M,,and M,, 


Myx—y 


» 


Fic. 1. 


and the torques, M,, = — M,, are given per unit length of side. 
We may also define the curvatures, a typical value of which is 


(dw/ dz) 
a1) 


Piss Ox? Ox 


Then (see footnote 2) we may show that the following tensors occur 
in the theory of moment and torque loaded thin plates with small 
deflections: 


moment tensor 


(12) 


M,, —M,, ) 


curvature tensor 


ow 1 dy) dz) 
ax? 2 ax ay 
ia | 1 dz) Q(aw/ ay) ew 
2 Oy ax 


(13) 


: 
M 
Shy J, 
Z Jy 
ax 
ae 
zz fy 
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wie r yz r 
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It will be noted that these are tensors, in the sense of Eq. 3, and 
are symmetrical. Also, because they are tensors they satisfy all the 
transformation, invariance and other relations governing tensors. In 
particular, because they are symmetrical they may be put in principal 
axis (that is, diagonal) form. Also, they satisfy the Mohr Circle 
transformation construction and the usual rotation relations for two- 
dimensional tensors. 

It may be shown that the equations governing the action of thin, 
laterally loaded small deflected plates are (see footnote 2) : 


equilibrium equations 


div 

( )( M.; 

ox ay 
ou NV = (M) 


—M., é ) 
M,, \M,, 


or 


boundary conditions 


in which M,, and M,, are the edge moments. 


Relation Between Curvatures and Moments (Hooke’s Law) 


It may be shown (see footnote 2) that the following relation holds 
between and 
Vv 
— —— I, 8s. (18 
DG — y*) 18) 


In Eq. 18 9, is the trace or first invariant of the moment tensor and 
D is the plate stiffness. 
Finally, there is a 


compatibility of curvature relation 


oy? Ox? Oxdy 
THE ANALOGOUS QUANTITIES AND RELATIONS FOR THE TWO-DIMENSIONAL THEORY 


In view of the foregoing, it follows that the following analogous 
quantities and relations may be set up for the two-dimensional isotropic 
elasticity-thin plate fields: 


: 
Git 
(14) 
z 

c 
iy 
(16) 
: 
or 
: 
1 
( / ) 
‘i 
ay 
P 
"hh 
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the tensors 


K is a constant in the above. Also 


n ~h@ 
(21) 


or 


Eq. 2 ~ h-Eq. 13] 


equilibrium equations 


div T = 0 ~ div m = 0 


Eq. 5 ~ Eq. 14 


boundary conditions 


TN =X ~ MN = M| 
or (23) 
Eq. 7 ~ Eq. 16 


Hooke's Law 


Eq. 9 ~ Eq. 18 


compatibility condition 


Eq. 10 ~ Eq. 19. 


We see therefore that there is a direct analogy between the two- 
dimensional equations of the linearized theory of isotropic two-dimen- 
sional elasticity and the governing equations for the thin, small deflec- 
tion theory of moment and torque loaded plates. 

Because of the analogy, it follows that solutions in the two-dimen- 
sional theory of isotropic elasticity may also be taken as analogous 
solutions in the theory of moment and torque loaded thin plates. In 
applying the analogy we compare a plate in the x, y plane (deflecting 
laterally in the s direction) to a two-dimensional deformable body 
loaded by stress or forces in the x, y directions. In particular, the 
analogous quantities are (using the subscript £ for elasticity solutions 
and P for plate solutions) : 


289 
T ~ 
K h? 
(20) 
or 
1 
= 
ie 
(22) 
or 
(24) | 
25 
(Zo ) 
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Elasticity Thin Plates 


~ 


My My, 


h 


Illustrative Example 


As an illustration of the application of the analogy, we consider the 
following problem in the theory of elasticity, Fig. 2.5 


| 

Cc 
i” 

= 


Fic. 2. 


°S. TIMOSHENKO AND J. N. Goopter, ‘“‘Theory of Elasticity,” second edition, New York, 
McGraw-Hill Book Co., Inc., 1951, p. 32. 


: 
: 
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The beam shown, with stresses 


ds(x*y — 
3d (27) 
— 


II 


oy 


Try 


satisfies all the requirements of the two-dimensional linearized theory of 
elasticity. Along the top and bottom sides the normal stress ¢, is 
uniformly distributed. On the side x = / the normal stresses follow a 
linear and cubic parabolic variation. The shear stress on the top and 
bottom is proportional to x and is parabolic on the face x = /. 

Thus, the solution shown represents a thin plate solution, with 


M,./ Kh’ in place of 
M,,/Kh? in place of ¢,! 


(28) 


M., 


in place of 7.,. (29) 


From Hooke’s Law the elasticity strains may be obtained as 


ez = 


&y = 


= 


and therefore the radii of curvature for the plate are known, using 
Eq. 26. 


THE THREE-DIMENSIONAL FORM OF THE ANALOGY 


We consider next a thin plate having small deflections and subjected 
to transverse and shear loads in addition to the edge moments and 


torques, see Fig. 3. 
Then, it may be shown that the equilibrium equations for this 


structure are 


, 

Ox oy 

oM,, , 
Ox oy 


, 
Ox oy J 


= 0 


(31) 
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which is also given by 


div m = 0 


aa 
Ox OY OZ 
—Q.2 


In the above a new tensor, 91, defined by 


has been introduced. This tensor represents the three-dimensional 
generalization of the two-dimensional form given in Eq.*12 and in- 
cludes the effect of transverse and shear loads. 


We are therefore led to introduce a generalized three-dimensional 


| 
3 Myx 
ox 
Onto 
Fic. 3, 
| 
| 
(32) 
or 
( -My —-Q.2 
 -0.: 
M,. My, Q,2 (34) 
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This becomes, in its three-dimensional 


form of ®, as given in Eq. 13. 
form, 


2 “oy Ox oy 


yt and ® as formulated above are tensors 


It should be noted that 
in the sense of Eq. 3. 
Because of the generalized nature of these tensors one is led to 
postulate a three-dimensional Hooke’s Law as a direct generalization of 
the two-dimensional case, Eq. 18. This is given directly by 


In this equation, ® and ® are given by Eqs. 35 and 34, 9, is the trace 
of Mand &; is the 3 X 3 unit matrix. 
The boundary or edge conditions on the moments, torque and shears 


= (3) (37) 


are 


(38) 

-Q2 |Qz 


In Eq. 38, /, m and n are the direction cosines of a normal to the 
edge of the plate at the central plane z = 0, hence n = 0. The right- 
hand side terms are the applied edge moment, torque and shear. 

Also, the elements of the curvature tensor, ® satisfy the following 
six compatibility conditions 


: 
5 
1 1 1 
1 1 1 
1 
ar zy 2 
; 
; 
i 
: 
‘ 
; 
or 
| 
| 
: 
: 


ley 
Oz Ox 


We now turn to the relations and equations of elasticity. In the three- 
dimensional theory of elasticity, the fundamental stress and strain 


tensors are, 


Try 


(40) 


and 


1 / dav ou ) 
Ox oy 
1faow ou ) 
2 \ ax OZ 
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The equilibrium equations in three-dimensional elasticity are 


div T = 0 (42) 


(43) 


between the strain and _ stress 


Also, Hooke’s Law, the relation 
tensor is 


1 
(44) 


In Eq. 44, 9, is the trace of the stress tensor, 6; is the 3 X 3 unit 


matrix. 
The boundary conditions on stresses is given by 


(45) 


(46) 


satisfy the following six compatibility 


and finally, the strains must 
conditions, 


, Oe, OY yz 


dz? 


0°e, 
ix? 


a 
oy’ Oxdy 


(47) 
Oxdy as\ ax oy 
He, ( _ ote ) 


Ox oy 
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By comparing the above, it becomes clear that an analogy may be 
drawn between the postulated thin plate relations, Eqs. 31-39 and the 
elasticity relations Eqs. 40-47. 

This means that solutions in plate theory may represent solutions 
in the three-dimensional elasticity theory and vice-versa. 

In Timoshenko’s book (see footnote 3) many different plate solutions 
arefgiven. Each of these corresponds to a particular solution in a 
three-dimensional stressed structure providing it satisfies all the 
relations, Eqs. 31-39. The analogous quantities are the following: 


Elasticity Plate 
Kh? 
My 
Kh? 


Kh? 


K oh 2 


The procedure may be summarized as follows: 


1. Knowing w(x, y)p determine M,,, M,,, M.,, Q: and Q, using the 
known plate relations. These are all functions of (x, y) only. 

2. Using the values from (1), from Hooke’s Law (Eq. 36) determine 
the 1/r terms. Note, these are now functions of (x, y, 2). 

3. Check the compatibility equations, Eq. 39. If these hold, then 
the analogous relations, Eq. 48, follow. 


As an example, it is well known that the exact solution for a uni- 
formly loaded elliptical plate with built-in edges is given by 


296 
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2 2 2 
w(x, y) (1 (49) 


in which c is a constant and a, } are the major and minor elliptic radii, 
respectively. Using the known plate relations we obtain from Eq. 36 
the various elements of ® which when substituted in Eq. 39, give 6 
equations. The third, fourth, fifth and sixth of these are satisfied 
identically. The first and second will be satisfied if 


a = 6 (that is, for a circle) 


(50) 


and if 


Then using Eq. 44 we can obtain the strains and from these the 
deformations u, v and w; or, these could be obtained using Eq. 48. 


CONCLUSION 


Two separate analogies between clasticity theory-thin plate theory 
were developed. 

In the first, two-dimensional elasticity solutions are shown to be 
analogous to thin plates with small deflections and torque and moment 
loaded only. 

In the second, three-dimensional elasticity solutions are shown to be 
analogous to thin plates with small deflections and laterally loaded as 
well as torque and moment loaded. 

Examples of both were given. 


APPENDIX 


Symbols Used in This Paper 
T = stress tensor 
= strain tensor 

oz = normal stress in x direction 

oy = normal stress in y direction 

o; = normal stress in z direction 

Tyr = Tzy = Shear stress in x—y direction 

shear stress in y-z direction 


ae Using these we find the analogous stresses are given by 
12M,: 
oz h? 
i 
h 
1 
: 

: 

: 

: 
\ 
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= shear stress in s—x direction 
éx = normal strain in x direction 


é, = normal strain in y direction 


é, = normal strain in 2 direction 
Yry = Yyr = shear strain in x-y direction 
= Yzx = shear strain in y—s direction 
Yex = Yr: = Shear strain in z—x direction 
u = deformation in x direction 
v = deformation in y direction 
w = deformation in z direction, deflection in the z direction of thin plate 


v, Y, 3 = coordinate axes 
h = plate thickness 
v = Poisson's ratio 
E = Modulus of elasticity in tension-compression 


9, = the trace, the first invariant of the tensor 


100 10 
&3, S&2 = unit matrix, | 010 hos ) 
001 


M,, = plate moment per unit length on a plane perpendicular to the x axis 
M,, = plate moment per unit length on a plane perpendicular to the y axis 
Mey = — M,z = plate torque in the x—y direction 
q = lateral load per unit area on the plate 
Q, = plate shear force, per unit length, on a plane perpendicular to x axis 
Q, = plate shear force, per unit length, on a plane perpendicular to y axis 


a = moment tensor, defined by Eqs. 12, 34 


@ = curvature tensor, defined by Eqs. 13, 35 


typical curvature element of the curvature tensor 


D = plate stiffness, Dd.) 


X = boundary stress in x direction 
y = boundary stress in y direction 
Z = boundary stress in z direction 
M,; = boundary moment in the i direction 
M;; = boundary torque in the i—j direction 
Q; = boundary shear 


K =a scale factor constant in the analogy relations between 4 


and @ f 
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ANALOGY OF NONLINEAR SYSTEMS TO CLASSICAL DYNAMICS * 


BY 


RUEY-WEN LIU! AND GILBERT H. FETT? 


ABSTRACT 


The two-dimensional autonomous system in the general form is considered. 
The solution is in terms of the phase trajectories in the phase plane. In this paper, 
a force field over the phase plane is formulated such that the path of a particle in this 
force field coincides with the path of the trajectories, provided the initial condition is 
properly chosen. Hence, the trajectories can be obtained by solving the path of the 
particle subject to the formulated force field. It is further developed that, for any 
given system, it is theoretically possible to formulate a conservative force field to give 
the same property. Consequently, it is found out that the trajectories can be solved 
either analytically by a simple substitution or graphically by means of the family of 
trajectories which are orthogonal to the trajectories of the original system. Various 
examples including the generalized Van der Pol equation are given gs illustrations. 


I. INTRODUCTION 


Equations under consideration are a pair of first order autonomous 
differential equations, defined as System (1), in the general form 


= p(x, y) (1) 
y = q(x, y), 


where the dot denotes the time derivative. It is assumed that the 
partial derivatives of p(x, y) and q(x, y) with respect to x and y exist 
and are continuous in the domain considered. If p(x, y) = y, then 
System (1) becomes a single differential equation of second order: 
= q(x, 

If the independent variable time ¢ is eliminated, System (1) can be 
written as 


dy _ q(x, 
(2) 


Since the partial derivatives of p(x, y) and q(x, y) exist and are con- 


* This paper is taken from a thesis submitted to the Graduate College of the University of 
I}linois in partial fulfillment of the requirements for the Ph.D. degree in Electrical Engineering 
by the first named author under the direction of the second. The assistance of a General 
Electric Fellowship is gratefully acknowledged. 

1 University of Notre Dame, Notre Dame, Ind. 

2 Department of Electrical Engineering, University of Illinois, Urbana, III. 

3 The boldface numbers in parentheses refer to the references appended to this paper. 
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solution of (2) exists with the form 


where c is an arbitrary constant. Equation 3 can be represented as a 
family of curves with parameter c in the x—y plane which is called the 
phase plane. The curves in the phase plane corresponding to (3) are 
called the phase trajectories, or simply the trajectories. Equation 2 is 
the differential equation of these phase trajectories. 

The purpose of this paper is to formulate a force field such that an 
analogy between 


(A) the trajectories of System (1) in the phase plane, and 

(B) the orbits of a particle with unit mass subject to the formulated 
force field are made. Hence the trajectories of System (1) can be 
obtained by solving the latter. 


The method of analogy is given in Section II]. After the analogy is 
made, several special classes of p(x, y) and q(x, y) are considered and 
solved in Sections III and IV, and the functions p(x, y) and q(x, y) in 
general are considered in Section V. Examples, such as the general- 
ized Van der Pol equation, # + f(x)# + g(x) = 0, are used as illus- 
trations. 

Il. THE ANALOGY 


When the time derivative of both sides of System (1) is taken, one 
obtains 

Op 

oy 


q = P(x, y) (4a) 


and 


4+ = O(x, (4b) 


oy 


The right-hand sides of (4a) and (40) are functions of x and y only, 
hence they can be defined as P(x, y) and Q(x, y). 
Let the force field F(x, y) be defined as 


F(x, y) = 4P(x, y) + 9Q(x, y) = + Jy. (5) 


By doing so, Eq. 5 becomes the equation of motion of a particle with 
unit mass subject to a force field whose x- and y-components are P(x, y) 
and Q(x, y). The force field so defined is called the associated force field 
of System (1), and Eq. 5 is called the associated equation of motion. To 
be distinguished from the trajectories, the paths of the unit mass particle 
are called the orbits. The associated momentum field M(x, y) of the unit 
mass particle is thus defined as 


F(x, y;c) = 0, (3) 
i 
: 
: 
th 
& 
‘i 
i 
7 
Ox 
Z 
= 
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(6) 


M(x, y) = ip(x, y) + jq(x, = + jy 


which is equivalent to System (1). 

Since the phase trajectories of System (1) should satisfy (6), and 
since (6) is a special integral of (5), the solution of System (1) is the 
subset of the general solution of (5). Furthermore, this subset should 
also satisfy (6) at all times, in particular, at time zero. Therefore, the 
following analogy can be stated: 


(A) The family of the phase trajectories of System (1) is a sub- 
set of (B) the totality of the orbits of a unit mass particle subject 
to the associated force field of System (1). Furthermore, the in- 
itial condition of this subset should satisfy the associated momen- 
tum field. 


It is noticed that, after the initial time, these trajectories also satisfy 
the associated momentum field. 

By this analogy, a purely mathematical problem has been changed 
toa dynamical one. Since the physical quantities have been attributed 
to facets of the mathematical problem, intuition may enlighten the 
solution to the problem. Furthermore, for certain types of force field 
obtained by (5), the solution of System (1) may be obtained in a straight- 
forward manner. 


III. A CLASS OF p AND gq 


Consider the class of p and g of System (1) of which the associated 
equation of motion is integrable to the form 


(7) 


f(z, x, = 0, 


where c is an arbitrary constant. Since (7) should be satisfied by all 
the orbits which are the solution of the associated equation of motion, 
it should also be satisfied by the trajectories. Hence the trajectories 
should satisfy (1) and (7) simultaneously. If (1) and (7) are inde- 
pendent, then « and y can be eliminated between them, and the phase 
trajectories of System (1) can be simply written out as 


flb(x, a(x, 956] = g(x, 0) = 0. 


Hence the following theorem is obtained : 


Theorem I: \f the associated equation of motion of System (1) can be 
integrated to the form of (7), and if (7) and System (1) are independent, 
then the phase trajectories of System (1) can be obtained by a simple 
substitution, 
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The following examples will illustrate this point. 


Central Force Fields (2) 


If the components of the force field P(x, y) and Q(x, y) satisfy the 
condition 


yP (x, y) = xQ(x, ¥), (8) 


then the 6-component of the force field vanishes identically. This force 
field is called the central force field. Consequently, the equation of 
motion can be integrated to give the equation of conservation of angular 
momentum, that is, in the rectangular coordinates, 


xy — ty =. (9) 


The associated force field of the following example is central. 
Example 1. Consider the system of undamped harmonic oscillator : 


Zz=y 
(10) 


where A is any positive constant. The components of the associated 
force field are obtained from (4): 


P(x,y) = — Ax 
Q(x, y) = —Ay 


which satisfy (8). Hence the associated force field is a central one. 
Consequently, the angular momentum of the unit mass particle is con- 
served. The trajectories of System (10) thus are obtained by substi- 
tuting (10) into (9): 

Ax? + y? 


For positive A, the trajectories are ellipses in the phase plane, which are 
closed trajectories; hence System (10) will oscillate in such a way that 
the particle moves on the constant angular momentum curve. 


Conservative Force Fields (3) of Type I 


Suppose the components of the force field P and Q satisfy the 
relation 


(11) 


which is equivalent to Curl F(x, y) = 0, then the force field is conserva- 
tive. Consequently, the equation of motion can be integrated to give 
the integral of energy: 


+ + Vix, y) = E, (12) 


PY 
: 
ely 
= 
: 
: 
: 
Bick 
4p : 
of 

dy Ox 


NONLINEAR SYSTEMS AND CLASsICAL DyNAMICS 303 


Oct., 1961.] 


where £ is the total energy. The potential energy V(x, y) is given by 


Vix, = — f P(x, yo)dx — f “O(x, ydy (13) 


where (Xp, Vo) is the arbitrarily chosen reference point. Furthermore, if 
Eqs. 12 and 1 are independent which is defined as Type J, then the 
trajectories of System (1) can be obtained by a simple substitution. 
The following example will be given as an illustration. 


Example 2: Consider the system 


t=y 


y=x+ Bel*+A — (14) 


where A, B and k are any constants. From (4), the following com- 
ponents of the associated force field are obtained, 


x + + A — (§)ky’ 
= y — Aky + (3)k*y? — key, 


P(x, y) 
Q(x, ¥) 
which satisfy (11). Hence the associated force field is conservative, 
and the trajectories should satisfy (12) as well as (14). If the origin 
is chosen as the reference point, the potential field is obtained from (13) : 


B 


Substitute this potential field and (14) into (12) to give 


(})A* + + Bx — ($)kBy* + (2) =E 


which is the exact solution of System (14) in the phase plane. Here the 
total energy E is the parameter of the family of the phase trajectories. 

It is noticed that these two examples are not solved by the con- 
ventional mathematical way. After the analogy is made, through the 
knowledge of classical dynamics, they are solved in a straightforward 
manner. There are still many other force fields which make the equa- 
tion of motion integrable, but they will not be considered here. 


IV. THE cass: = 29 
oy Ox 


Conservative Force Fields of Type II 


The case for which the associated force field is conservative and the 
integral of energy (12) and System (1) are dependent is defined as 
Type II. ‘The following theorem is thus introduced, 
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Theorem II: The necessary and sufficient condition for which the asso- 
ciated force field is a conservative field of Type IT is 


Op 0g (15) 
ay 


Proof: If (15) holds, then, from 


proved : 


(4), the following relation can be 


oP 


Oy Ox” 


Hence the associated force field is conservative. The potential energy 


V(x, y) is obtained from (13) : 


V(x, y) = — 3[p2(x, + @ix, + Vo, 


where 


Vo p? (xo, Vo) 4 (Xo, yo) | 


and (xX, Yo) is any arbitrarily chosen reference point. The quantity Vo 
is defined as reference potential which is constant for any chosen reference 
point. If (16) is substituted into (12), the integral of energy becomes 


(18) 


The necessary condition for which the solution of these simultaneous 
Eqs. 1 and 18 exists is 


E = Vo. (19) 


Under this condition, System (1) and the integral of energy (18) are 
dependent equations. 
The reverse part is proved by use of the relationships 


(20) 


Since Eqs. 1 and 12 are dependent, when # and y are eliminated be- 
tween them, V(x, y) has to satisfy the equation 


(P+¢@)+ Vix,7) =E 


identically. If (20) is applied, the same result is obtained from both 
of these equations, namely, condition (15). This completes the proof. 
It is noted that Eq. 19 gives the condition that the particle moves 


: 
4 
: 
| (16) 
: 
3 
P(x,y) = —— V(x, y) 
Ox 
O(x,y) = —— Vix,¥). 
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with the same constant energy, that is, Vo, along any trajectory in the 
domain where (15) is satisfied. This is the case because the initial 
condition is pre-determined by System (1). 

Since the integral of energy and System (1) are dependent equations, 
the trajectories can not be obtained by a simple substitution as in the 
case of Type I. Other methods should be further developed. One of 
them has been completed and is given in the following paragraphs. 


The P-Traj 
The family of trajectories which are orthogonal to the family of 

phase trajectories of System (1) is called the perpendicular trajectories, 

abbreviated as ‘‘P-traj.’’ In contrast to the P-traj, the trajectories of 

System (1) are called the actual trajectories. 

The differential equation of P-traj is given as 


_ P(x, 9) 
q(x, 


which is the same type of equation as (2). Hence the solution of (21) 
exists. Furthermore, if p and q satisfy (15), P-traj can be obtained by 
a simple integration. 


Theorem III: Assume that condition (15), ap/dy = dq/dx, is satisfied. 
Then the P-traj is given by 


U(x, y) = yo)dx + [a y)dy (22) 
70 yo 


where c is the parameter of the P-traj. 


Proof: Equation 15 is equivalent to Curl M = 0. Hence there exists 
a scalar function U(x, y), defined by (22), such that (3, pp. 279-280) 
Grad U(x, y) = M(x, y). Thus the direction of M is perpendicular 
to the constant U-curve. Since the direction of M is the same direction 
as the actual trajectories, the actual trajectories are also perpendicular 
to the constant U-curve. This completes the proof. 


After the P-traj is obtained, the actual trajectories can be obtained 
by tracing a path normal to the P-traj. Hence the actual trajectories 
are solved graphically. 

Example 3: Consider the system 


&¢ = Ay + g(x) 
Ax + 


| 
| 
| 
| 
—_— 
ax | 
| 
| 
| 
| 
| 
| 
| 
7 
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where A is any constant. This system may represent a two-loop elec- 
trical circuit coupled by a mutual resistance. Both loops have identical 
self-inductances but different arbitrary nonlinear resistances as shown 
in Fig. 1. These nonlinear resistances could be tubes, transistors or 


NL 2 NL 
RESISTOR A RESISTOR B 


Fic. 1. Two-mesh network with nonlinear resistors to illustrate Eq. 23. 


other devices. The variables x and y here are the two loop currents 
7, and 722. 

It is seen that System (23) satisfies (15). Hence, the associated 
force field is a conservative field of Type II]. The P-traj is obtained 
from (22): 


U(x, vy) = A(xy — + seoax + h(y)dy =c (24) 


which can be plotted if g and / are given. 
Consider the simple case: A = 1, g(x 0, and h(y) = (3 
System (23) becomes 


#— —x = 0. 
The P-traj of (25), for (xo, yo) = (0, 0), is thus 


( 
U(x, y) =xy+-— 
6 


which is plotted with the actual trajectories in Fig. 2. 

It is noticed, from the previous example, that even though the 
original system is not conservative, the associated force field may be 
conservative. 

It can be concluded that if the associated force field is conservative, 
then the trajectories can be obtained either (A) analytically by a 
simple substitution in the case of Type I, or (B) graphically by means of 
P-traj in the case of Type II. 


9 
(25) 
: 
: 
: 
eed 
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V. THE CONVERSION OF NONCONSERVATIVE FIELDS TO CONSERVATIVE FIELDS 


Converting Factors 
System (1) in general is considered in this section. 
Multiplication of the right side of System (1) by 7 (x, y) gives the 


new system 


q(x, y)T (x, 


Fic. 2. The actual trajectories (dotted lines) and the perpendicular trajectories (solid lines) 
of the equation ¢ — 442 — x = 0 corresponding to the network of Fig. 1. 


The differential equation of trajectories of the new system is 


dy _ y) 
dx p(x, y) 
which is the same as that of the original System (1). These two sys- 


tems give the same trajectories. Hence, as far as trajectories are con- 
cerned, System (1) can be replaced by the new System (26). 


307 
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If 7(x, y) is chosen such that the relation 


ay p(x, T(x, ¥) = ax y)T (x, y) 


holds, then System (1) is converted to a new System (26) of which the 
associated force field is a conservative one of Type II. The function 
T(x, y) which satisfies (27) is called the converting factor. 

Since (27) is the partial differential equation of integrating factors 
(4) of the differential equation, g(x, y)dx + p(x, y)dy = 0, and, since 
there exist infinitely many integrating factors for this type of differen- 
tial equation, there exist infinitely many converting factors. 


Theorem IV: There exist infinitely many converting factors 7(x, y). 
Consequently, the original System (1), theoretically, can always be 
converted to the new system 


& = p(x, y)T (x, y) 
y = q(x, y)T(x, y) 


of which the associated force field is a conservative one of Type II. 
Hence the actual trajectories can be obtained graphically. 


REMARK: The multiplication by y) may introduce 
tinuities, at points (x,y) which are discontinuities of 7 (x,y) or 
extraneous solutions (3, p. 441): curves along which Te) = = (), 
These complications are not considered here. 

Now, the problem of converting a given system to a new system, of 
which the associated force field is conservative, depends upon finding 
one of the infinitely many converting factors. This is essentially a 
problem of finding any solution which satisfies the partial differential 
equation (27). This is theoretically possible, but, generally speaking, 
it is rather difficult. However, for certain p and gq, the converting 
factor can be obtained readily. This is illustrated by the following 
examples. 


Example 4: Consider the generalized Van der Pol equation 
f(y)y + = 0. (28) 


Let F(y) = f f(y)dy and x = y + F(y). System (28) can be trans- 


formed to the following expression : 

— gly) 
(29) 
y=x — F(y). 


This x-y plane is called the Liénard plane (s). 


J. 
: 
(27) 
3 
i 
: 
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The converting factor of System (29) can be obtained by solving 
(27). Assume that 7(x, y) is a function of y only, then (27) can be 
simplified to the form 


d 
dy = T(y) 


which can be solved to give 


G(y) 
g(y) 


T(y) = (30) 


where G(y) = Exp (- f - dy). Consequently, the associated force 


field of the new system 


is a conservative field of Type II. The P-traj is obtained from ( 


— xG(y) — J [ro Cy) 


g(y) 


Jay =h (32) 


where h is a constant which depends on x) and yo. Equation 32 is the 
P-traj of the generalized Van der Pol equation in the Liénard plane. 
For given f and g, P-traj can be plotted, and hence the trajectories can 
be solved graphically. 


Example 5: Consider the system 


F(t) +x =0. (33) 


This equation could be an electrical circuit consisting of a nonlinear 
resistance in series with a linear inductance and a linear capacitance, 
or it could be a mechanical system with a nonlinear damping and a 
linear spring force. 

It is seen that (33) is a special case of System (29) with g(y) = y. 
Hence the converting factor can be obtained from (30): 7(y) = 1/9”. 
The P-traj is obtained from (32) : 


( 
| 
| 
| 
| 
| 
| 
| 
G(y) (31) 
= [x — F(y)]— 
| 
| 
| 
ie 
| 
| 
| 
| 
| 
| 
vf F(y) (34) 
dy = h. : 
| 
| 
2 | 
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The trajectories of System (32) thus can be obtained graphically by the 


use of the P-traj. In the following example, the procedure has been 
carried out. 


The Van der Pol equation. Consider the Van der Pol equation 
=0, (35) 


which is a special case of (28) with f(y) = — a(1 — y*) and g(y) = y. 


13 
The P-traj is found from (32) with F(y) = — a (> _ * : 


2 


x 


Fic. 3. The actual trajectories (dotted lines) and the perpendicular trajectories (solid lines) 
of the Van der Pol equation (35) with a = 1. 


‘ 

: 

: 

| ) =k (36) 

£ 

x 

4 

3 

: 

-3 -2 2 
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where k is a constant which includes x9 and yo. This is the P-traj of the 
Van der Pol equation in the Liénard plane. For a = 1, the P-traj and 
actual trajectories are plotted in Fig. 3. 

It is seen from Fig. 3 that every actual trajectory eventually ap- 
proaches a closed trajectory. This closed trajectory is called a limit 
cycle. The limit cycle corresponds to a steady state oscillatory solu- 
tion in the actual system. Other trajectories which approach the limit 
cycle are the transient solutions. 

There are two points which are noticed in Fig. 3 and which may be 
further investigated : 


1. It seems that all the trajectories outside the limit cycle approach 
the maximum curvature curve of the P-traj first. Then they are led to 
the limit cycle along this maximum curvature curve. This fact may 
indicate some information about the transient behavior of this system. 

2. It is seen from (36) that the P-traj varies continuously with the 
parameter a. Hence, with a varying, the transition of the P-traj can 
be readily found. The transition of the actual trajectories thus can be 
studied through the orthogonal relationship. 


VI. CONCLUSIONS 


Summary 


An analogy from a mathematical problem to a problem of classical 
dynamics is developed in this paper. 
After the analogy is made, three results are given: 


1. If the associated equation of motion of System (1) admits an 
integral in the form 


f(z, y;c) = 0 


and if the integral is independent of System (1), then the phase trajec- 
tories of System (1) can be obtained by a simple substitution. Un- 
fortunately, only a few p and q belong to this class. 

2. If the relation (15), ap/day = dq/dx, is satisfied, then (a) the 
associated force field is conservative and (b) the trajectories can be 
obtained graphically by use of P-traj. 

3. If relation (15) is not satisfied, that is, considering p and gq in 
general, then, theoretically, there exists a converting factor which con- 
verts System (1) to a new system which satisfies the relation (15). 
Hence results given above again will be obtained. The converting 
factor of the generalized Van der Pol equation has been found. The 
P-traj and the actual trajectories of the Van der Pol equation are 
plotted and discussed. 
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Advantages of the P-traj Method 


There are three advantageous points about the method of P-traj. 


1. If the P-traj is found, then the trajectories can be obtained 
graphically in an easy and accurate manner. 

2. Because of the orthogonal relationship between the P-traj and 
the actual trajectories, the qualitative behavior of the latter could be 
analyzed by studying the curvatures and the relative positions of the 
former as suggested in the example of the Van der Pol equation. 

3. It seems possible that the problem of finding the proper functions 
f,"2, h or F of Systems (23), (28) and (33) to give prescribed trajec- 
tories can be solved by means of the P-traj._ If the actual trajectories 
are given, the required P-traj can be obtained through the orthogonal 
relationship. It is seen from (24), (32) and (34) that the P-traj are 
explicit functions of f, g, # and F. Hence the most proper functions 
f, g,h and F could be found to give the P-traj which is closest to the 
required P-traj. 


Further Suggestion 


It is noticed from Theorem IV that every System (1) can be con- 
verted to a new system of which the force field is conservative. How- 
ever, the properties of the conservative field have not been used exten- 
sively in this paper. It might be more profitable if the knowledge 
concerning the conservative field, such as the Hamiltonian Principle, 
the contact-transformations, the least action principle, etc., is further 
used to solve the problem. Conversely, some knowledge about the 
trajectories in the phase plane may lead to the knowledge about the 
orbits of a particle in a conservative field. 
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READER OPINION 


ADDENDUM TO “PONTRIAGIN’S MAXIMUM PRINCIPLE AND 
THE PRINCIPLE OF OPTIMALITY” * 


BY 
C. A. DESOER'! 


The purpose of this note is to show that Eqs. 9 of the paper reduce 
in all cases to the form 


J 


( 
dt j=0 Ox; (1 ) 


It has been shown in the paper that if the optimal control u(x) is 
an interior point of the set 2, then Eqs. 9 reduce to Eq. 10. Consider 
now the case where for all states in some small neighborhood N(x) of 
the state x, the optimal control u° is a boundary point of 2. For con- 
venience, let 


$(x, u) = (g(x), f(x, u)) = ¥ g;(x)f;(x, u). 
7=0 


Since u®(x) is a boundary point of Q and since u°(x) is the point at 
which (x, u) takes its maximum value (over the set Q), it follows that 
V.ud(x, u), the gradient of ¢(x, u) with respect to u, is parallel to the 
outward normal of 2 at u°(x). Consider now the point x + dx; where 
dx; is parallel to the 7‘ axis of coordinates and sufficiently small so that 
x + dx;eN(x). Consequently, u°(x + dx;) is on the boundary of Q 
and u°(x + dx,;) — u°(x)«<du°/dx,; is parallel to the plane tangent to 
the boundary of 2 at u°(x). Therefore, the vector du°/dx; is orthogonal 
to the normal at u°(x) and to V.@ evaluated at (x, u°(x)); hence, its 
scalar product is zero, viz. 


n af; 01 0 


t=1 \ J OX; 


Therefore, Eqs. 9 reduce to Eq. 10. 

This geometric reasoning is an interpretation of a similar conclusion 
reached through variational techniques by Dr. L. Berkovitz in a Rand 
Corporation report P-2306 dated May 12, 1961. 


* Jour. FRANKLIN INsT., Vol. 271, pp. 361-367, (1961). 
1 Department of Electrical Engineering, University of California, Berkeley, Calif. 
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THE FRANKLIN INSTITUTE 


INSTITUTE ESTABLISHES SCIENCE INFORMATION CENTER 


The Library of The Franklin Institute has inaugurated a special 
Science Information Service (SIS) designed to aid industry and scien- 
tific institutions in research and development projects. The chief 
purpose of SIS is to aid industrial technical staffs and organizations in 
planning research, by reducing laboratory time and effort and by avoid- 
ing duplication in research effort. 

SIS is essentially a special library research service available to 
government, business, industry, and individuals who need scientific 
literature searches done by experts in the field, and who recognize the 
fact that such a service may save firms and research staffs hundreds 
of hours of individual reading time and perhaps thousands of dollars 
in avoiding laboratory duplication. 

The Service, which offers literature searches in any area in mathe- 
matics, engineering, physics, chemistry, and industrial processing, 
permits organizations and individuals to keep abreast of developments 
in their respective fields. 

SIS provides literature searches on specific scientific and technical 
subjects to the extent indicated and prescribed by inquirers. Sub- 
divided, the searches may include: (1) bibliographic compilations 
(author, title, source); (2) annotated bibliographies; (3) abstracts of 
articles on a regular basis; (4) current awareness, to keep research 
staffs advised of the most recent developments in their scientific inter- 
ests; and (5) translations of foreign technical reports and journals. 

The abstracting of articles ((3) above) is designed to cover exhaus- 
tively a specific field where current scientific information is required, 
and the service is carried out through scanning of available world 
literature related to the client’s topic with a view to providing 250 to 
500-word abstracts. 

Current Library resources used by SIS include more than 200,000 
reference works on all phases of engineering, physics, chemistry and 
mathematics, some 5,000,000 U.S., British and Swiss patents, plus 
national and foreign indices and abstracting services. Resources also 
include more than 2,000 currently received scientific journals, about 30 
per cent of which are from countries abroad, and a sizable collection 
of technical reports. 

Readers interested in SIS may write for further details to the 
Science Information Service Head, The Franklin Institute, 20th and the 
Parkway, Philadelphia 3, Pa., U.S.A. 
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THE FRANKLIN INSTITUTE LABORATORIES 
J. G. RICHARD HECKSCHER, ACTING DIRECTOR 


MACROSCOPIC FERROMAGNETIC RESONANCE EQUATION 
FOR ARBITRARY MAGNETIC FIELD * 


BY 


C. J. LIN' 


ABSTRACT 


Most ferromagnetic resonance experiments are arranged in such a way that the 
Kittel resonance equation can be applied in its simplest form. The physical implica- 
tion of this requires that: (i) the externally applied static magnetic field Ha be large 
enough that the magnetization M of the sample be parallel to the applied field; (ii) 
the applied field and the sample magnetization be along one of the principal direc- 
tions of the anisotropy energy; and (iii) smail angle precessions of M about the static 


equilibrium direction be assumed. Occasionally, cases are encountered which do not 
allow such a simplification. The present work indicates that the equation for the 
more general case, still assuming small angle precessions, can be derived in a manage- 
able form by an orthogonal transformation in real space. The various energy terms 
involved in the transformed equation are then treated separately. 


I. INTRODUCTION 


A resonance equation describes the favorable condition for resonance 
absorption. Unless higher mode spin-wave excitation or high damping 
phenomena are considered, small oscillations and negligible damping 
torques are invariably assumed in the derivation of a resonance equation 
for the fundamental mode. The usual procedure assumes the mag- 
netization vector M to have sinusoidal components M;,’ exp (wt) and 
M,/ exp (twt — ir/2) along the mutually perpendicular directions & and 
n, respectively, these directions being perpendicular to the static equi- 
librium direction @, (Fig. 1). Substitution of these components in the 
equation 


M = yM xX H (1) 


vields a pair of linear homogeneous equations in M;’ and M,’, where y 
is the magnetomechanical ratio and H is the sum of the externally 
applied static magnetic field H,, the demagnetization field H,z and the 
anisotropy field H,. From the secular equation, one obtains the pre- 
cessional frequency w in terms of H, and other parameters. ‘This is the 
original Kittel (1)? equation for complete line-up (that is, parallelism 
of M, H,, Hu, H. is implied). 

For less restricted conditions, the resulting pair of equations may not 
be homogeneous in M,’ and M,’. Under these circumstances, the 


* This work has been supported by Wright Air Development Division under Contract 
No. AF (616)-7499. 
' Magnetics and Semiconductors Laboratory. 
* The boldface numbers in parentheses refer to the references appended to this note. 
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following alternative set of equations of motion can be used, 


” 


bo 


where £ is the total free energy and ¢:, ¢, are the projections of the 
precessional angle on the ¢é and nf planes, respectively. A Taylor ex- 
pansion of E in terms of ¢:, ¢, around the static equilibrium direction 
yields a pair of first order homogeneous equations of which the secular 


equation is (2) 


Fic. 1. 


In using this expansion and neglecting higher terms, it is tacitly 
assumed that the amplitudes of ¢; and ¢, are approximately the same; 
hence, it is not, in general, correct to wite 


M?sin 6l 06 d¢? |’ 
unless the pole of the spherical polar coordinates is nearly perpendicular 
to M. 

Extension of the Kittel equation of all directions of M within 100 or 
110 planes of cubic crystalline anisotropy was done by Bickford (3). 
More general expressions for arbitrary M direction were derived 
quantum mechanically by Van Vleck, using a quadrupole-quadrupole 
term in the Hamiltonian (4). They were also derived classically by 
MacDonald, using a variation method, for M, H, line-up condition (5) ; 
that is, H, is assumed to be so large that M is practically parallel to H.,,. 


: 
: 
/ 
\ 
\ 7a 
J 

: a 
: 2 
dn 

al 


Oct., 1961.] THe FRANKLIN INstITUTE LABORATORIES 317 


For most experiments, 3-cm. band frequencies are used. The H, for 
free electron spin resonance is approximately 3000 oersted, while H, 
can be of an order of 10° if M ~ 10? and 2K ~ 10°. Hence, the angle 
between M and H, may be considerable and M may not even be able 
to point close to the vicinity of hard directions. 

The following discussion is a way of handling the more general 
situation without imposing severe limitations or applying drastic 
assumptions. 


Il. MORE GENERAL APPROACH 
Let @ be the unit vector along M and a; its components or the direc- 
tion cosines. To find the static equilibrium direction ¢ of a, it may 
be stated that 


bE = = 0 (5) 


subject to the constraint 


2= 2a-ba = 0). 


— = 0. (7) 


Thus 


= hag. (8) 


of cubic crystal- 


For example, on differentiating the energy in the case 
line anisotropy, Eq. 8 yields 


MH,B; + M*Nia; — 2(K, + + = — da; (9) 


where 8; are the direction cosines of H, and N; the demagnetization 


factors. 
Equation 9 multiplied by a; and summed over @ is 


MH, a8: + M* Nia: — + Ke) Das 


+2K,Da¢=—X (10) 


where 8; are the direction cosines of H,. By observing the structure of 
Eqs. 9 and 10, a; can be solved for given 8; either exactly or more con- 
veniently by a few cycles of iteration, that is, the static problem can be 
considered as a separate one whose numerical solution can be found 


as accurately as desired. 
The principal crystalline axes x, y, z, are generally used as the coordi- 
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nate system for expressing the anisotropy energy, and consequently 
the total energy E of the system will be referred to these axes. The 
essential step is to transform this sytem into a primed one x’, y’, 2’, 
which is the same as &, 7, ¢. This can be done by two consecutive 
rotations, one bringing z into ¢ and the other x, y into &, ». The orthog- 
onal transformation matrix can easily be written as 


Cao + Saja3; — Ca, + Sarva; 


— sv(1 — a;?) 


V(1 — a;?) \(1 — a;*) 


— Sao + Caja; Sa, + (11) 
= — cv (1 a;") 
V(1 — 


V(1 — 


aA» a3 


where c and s are the cosine and sine of an angle yet to be determined. 


The direction cosines of M in the primed system are 


a = sin Q: 


a,’ = sin ¢, (12) 
[1 — (sin? + sin*® ¢,) 


The first and second derivatives of a;’ with respect to ¢;, ¢, are easily 

evaluated, and at the static equilibrium direction all those derivatives 
0a,’ 

are zero except — = =—— — = ——- = — 1]. Hence, ap- 
99, 

plying a; as the intermediate variables for the differentiation, Eq. 3 

takes the new form as the resonance equation: 


1 ak dF )( dE ) ( ek (13) 


Since — and y are presumably the principal directions of oscillation we 
must have 


= 1 and 


CE 


= (), (14) 


: 
Qe 
» 9 
¢ 0; — a3?) V(1 — a;*) 
¢ { ij Cc Q — 
— — a;*) 
| 
; 
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When the terms of (14) are written according to the transformation 
(11), there results an equation which serves to determine s (and hence 
also c = V1 — s?), In MacDonald’s derivation the internal fields 
obtained by a variation method were put into the form of effective 
demagnetization fields. Then, by dropping the off-diagonal elements of 
the demagnetization tensor, the resonance equation was obtained in the 
Kittel form. The present transformation has a similar effect, but it 
is achieved exactly and more directly. 

By means of (11), calculate the right-hand side of (13), term by 
term, for the standard forms of the free energy (6). 


(1) Magnetostatic Energy, ,.E 
= — MH, = — MH, ¥ (15) 


where 6; are the direction cosines of H,. Let 2”; ;, be the abbreviation 


The 


————. and E’ for 
0a;'da da 


E'y = — = — ad. (16) 


for 


3 


(17) 


(2) Demagnetization energy aE 


— aZN; (18) 


where NV; are the demagnetization factors (assuming the principal axes 
of the ellipsoid to be the same as the crystalline axes). Then 


Me A32N; M aZN; (19) 


yy = — A,2N; (cas + Saya)? ; 


+ (— ca; + Sa2a3)*N a3”)?N; |/(1 a3”) (20) 


= — Ao?N; Sa2 + Caya;)*N, 


+ (sa; + casa;)?N2 — c2(1 — as?)?N3]/(1 — as?) (21) 


ey = M AiiA oilV ; = M?[sc(N; N;) (a2? 


+ (s? — — Ni) ]/(1 — as?). (22) 
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(3) Cubic Crystalline Anisotropy 


where 


— (K,+ K;)/2 
Ls = Ks 3. 


L, Saand = Le be calculated separately. 


1 


oils’ = 4L, = 4L, > 


12Z, > 
12L,[ + — 

+ + + a2‘) |/(1 — a?) 
A 3i°A 2," 
12L,[ — — 

+ + + — a?) 


12L,[ 2sc Pas? — — (ai? — 


+ 2sca;?(ay4 + + ae) ]/(1 — a?) 


+ — a*) + 
30L» 
30L2[ — (ay? — + cas? 


= 30L.[ — 2" — — as”) scas? |. 


(30) 


(31) 


The less frequently encountered energy terms, such as crystalline 
anisotropy of uniaxial and hexagonal symmetries and magnetostriction 


energy, can also be treated in the same manner. 


A 
(23) 
Ki Ks 
: = + 
2 6 
(24) 
(25) 
: 
(20) 
E” 
= 
(28) 
29 
= 
= = 
2 


Oct., 1961.] TuHeE FRANKLIN INSTITUTE LABORATORIES 


III. APPLICATION TO SIMPLER CASES 


Case 1 


Consider an ellipsoidal sample with M and H, in the [001] plane 
of cubic crystalline anisotropy. From Eqs. 17, 22, 28, and 32 it is 
obvious that (14) is satisfied if s = 1,c = 0, ors =0,c¢ =1. Taking 
s = 1, then by Eqs. 14 to 32, the resonance equation (13) becomes 


w- 


7 (MH, a8: + + 


9 


4Li(a;4 = 6L2(a;® -+- MN; 4. | 
x (MH, > + M?(a,7N, >) 


— + a2‘) — 6L2(a;° + 
+ M?(a.N, aN.) + 24L +- 30L |. (33) 


Let a: = cos 6, a2 = cos 0, a; = 0; then Eq. 33 reduces to exactly 
the result obtained by Bickford in the M, H, aligned case. 


Case 2 


Spherical sample, H, and M in 110 plane of cubic crystalline anisot- 
ropy s = 1 also is a solution of Eq. 14, but not if these samples are 


1 
not spherical. Here a; = cos 6, a; = a2 = Va V1 — cos? 6. 
\ 
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x [ — 4L; — 6L2 a6 
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1 — a;? 


1 
M 
— 6L.(cos* 6 + } sin® @) + 18Z, 6 cos? 6 + 30L:2 cos? 
[M-H, — 4L;(cos‘ @ + 3 sin‘ 6) 

— 6L.(cos* + sin®@) + 6L, sin? — Lz sin*6]. (34) 


[M-H, — 4L;(cos‘@ + 3 sin‘ @) 


The crystalline anisotropy terms of Eq. 34 can again be reduced to 
exactly the same as Bickford’s. However, the terms (V,; — N;)M? and 
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(N. — N;)M?* cannot simply be added in Eq. 34 in order to include 
demagnetization energy for the general ellipsoidal sample since a more 
complicated solution than s = 1 is required to make gf"), + .E";'2' = 0. 


Case 3 
Spherical sample, first order cubic crystalline anisotropy, general M 


direction. The resonance equation (Eq. 13) for this case becomes 


Solving = 0, vields 


= sin] (36) 


9 


— a3?(ay? + + 


— | 


and c, which equals the cosine of the same bracket. Hence .,/’», 
ee” yy, edt’”’e» in Eq. 35 are given by Eq. 36 and Eqs. 25, 26, 27, re- 
spectively. 

IV. DISCUSSION 

In Sections IT and III it was shown that cases more general than the 
MH, line-up one can be treated without much difficulty. Solving for 
the static equilibrium first greatly simplifies the procedure. In the 
MacDonald or Van Vleck theory, the direction of the microwave mag- 
netic field was conveniently taken as the & axis (present notation). It’s 
magnitude was taken as negligible and did not enter the theory in any 
way. Hence the & direction seemed to be somewhat arbitrary. In the 
present work £ is fixed by Eq. 14. Actually, Eq. 13 alone will give the 
unique solution for the stable precession. 

Quantum mechanical justification of the macroscopic result can be 
made by following every step of Van Vleck’s derivation, except that the 
direction & is taken to be the quantization direction which is in general 
not the same as that of H,. After simplification of terms by various 
approximations as used by Van Vleck, there results 


+ 
boS; = as, + Ce (37) 


= w(— cS; — 


where a, },, be, C1, C2, and w are functions of direction cosines. Since 
{is presumably the direction of quantization, S; = 0;hencec, = c. = 0, 
By setting a = 0, the direction of & (and ) is fixed just as in Section II. 
hence the transformation relation is exactly the same as Eq. 11. Con- 
sequently a similar result can be expected. 
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YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 138-year history, there is vital need for 
The Franklin Institute effectively to promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive. They begin with students in 
the early grades of our elementary schools and continue throughout an individual’s 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased and increasingly effective. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 

If you wish to give securities: (1) Send the stock certificates to The Franklin 
Institute, Philadelphia 3, Pennsylvania, and leave the assignment space blank. (2) 
Execute a stock power and mail this to the Institute separately, leaving the assign- 
ment space blank. (The Institute, upon request, will gladly furnish blank stock 
power forms.) 
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NATIONAL BUREAU OF STANDARDS NEWS 


THE BLACK VOID REACTOR CONCEPT 


The National Bureau of Standards has been investigating a promis- 
ing concept in research reactor design—that of the Black Void Reactor’, 
or BVR. The BVR can simultaneously provide both fast and slow 
neutron fluxes that are spatially separated. An arrangement that 
provides this separation, conceived by C. O. Muehlhause of the Bureau 
staff, would permit radiation effects studies and thermal neutron beam 
studies (for example, neutron diffraction) to be conducted independently 
in the same reactor. In addition, this arrangement constitutes a type 
of flux trap (inverted) which yields higher fluxes for a given power 
input than does the usual research reactor. 

During a design study for the National Bureau of Standards Reactor 
(NBSR)?, to be constructed at the new NBS site at Gaithersburg, 
Maryland, consideration was given to the idea of incorporating a con- 
verter assembly into the NBSR design. This consideration led to the 
BVR concept of reactor design. Although the BVR appears to possess 
a number of advantages, it requires significant fuel element develop- 
ment and testing before it can be applied. Two new reactor types? 
under consideration by the Atomic Energy Commission may contribute 
technique to the development of the heavily loaded annular fuel 
elements required for the BVR. 

The general type of Black Void Reactor is the BVR-N, where N 
represents the number of heavily loaded (appearing black to thermal 
neutrons) cylindrical fuel assemblies in the configuration. As in other 
reactors, the fast neutrons released by fission in the fuel walls are slowed 
by the moderator, and a certain portion of these now slow, or thermal, 
neutrons return to the fuel assembly. Some of these thermal neutrons 
will be absorbed by the fuel walls and converted into fast fission neutrons 
to maintain the chain reaction. 

In the BVR, a portion of the fast fission neutrons and some neutrons 
of intermediate energy that have been partially moderated irradiate 
the void interior of the tube. A sizeable space is thus provided in 
which instrumented samples may be placed for radiation effects studies. 


1 For further information, see ‘“The Black Void Reactor Concept,” by C. O. Muehlhause, 
Nuclear Sci. and Eng., Vol. 7, p. 505 (June, 1960). 

2““High-Flux Research Reactor for New NBS Site,"” NBS Tech. News Bull., Vol. 45, 
(August, 1961). 

* High Flux Beam Reactor at Brookhaven National Laboratory, and High Flux Isotope 
Reactor at Oak Ridge National Laboratory. 
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At the same time, beam tubes in the moderator outside the fuel as- 
semblies can provide thermal neutrons for use in neutron diffraction 
studies. 

The simplest approach to the Black Void Reactor is the BVR-1, 
consisting of a single long cylindrical fuel assembly immersed in a large 
tank of moderating material. A practical example would be an as- 
sembly eight feet in length and eight inches in diameter. According 
to the design study, fast neutrons leaving the fuel walls would slow 
down in the moderator and return to the walls with a probability of 
0.543. The over-all multiplication of the system (number of fast 
fission neutrons produced per neutron absorbed in the fuel walls) must 
be adjusted to have the inverse of this probability value. At a power 
of three megawatts, a fast neutron flux of approximately 10'4/cm? sec 
would be available within the cylinder. A limited number of beam 
tubes may be placed in the moderator of the BVR-1. 

Another example of the BVR is the symmetrical BVR-3, which 
comprises three 6-in.-diameter, 8-ft.-long cylindrical assemblies spaced 
1.73 ft. apart. Here the return probability for neutrons leaving one 
assembly and being absorbed in any of the three assemblies is 0.642. 
Thus, the walls of each assembly must have a multiplication in excess 
of 1.506. 

The NBS design study showed that as a research reactor, the BVR 
should have the advantages of long assembly life (greater than 1 yr. 
at several megawatts) and safe operation due to long mean neutron 
lifetime. There is adequate excess reactivity for the thermal neutron 
beam tubes. The beams may be readily aimed at the space between 
assemblies and thus avoid virgin fission neutrons and gamma rays. 


TWO RADIOACTIVITY STANDARD MATERIALS 


Two additional standard samples of radionuclides—iron 55 and 
promethium 147—are being issued by the National Bureau of Standards. 
These two standards,‘ part of a group of 51 now available, result from 
the Bureau’s continual program to establish new and improved radio- 
activity standards. 

Both of the standards have wide applications for medical and indus- 
trial laboratory purposes. Iron 55, an electron-capturing nuclide with 
a half-life of 2.7 years, decays by X-ray and Auger-electron emission 
directly to the ground state. It is used extensively as a tracer both in 
biological and industrial research. 

Promethium 147 is a pure beta emitter utilized for generation of 


4Standard materials issued by the National Bureau of Standards are described in NBS 
Circular 552, Standard Materials. This Circular may be ordered for 35 cents from the Super- 
intendent of Documents, U. S. Government Printing Office, Washington 25, D.C. Up-to-date 
supplementary inserts, which are issued periodically, are available upon request directly from 
the National Bureau of Standards. 
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low-energy X-rays and energy calibration standards for X- and gamma- 
ray spectrometers. It has a half-life of 2.5 years. In addition to other 
research applications, this radionuclide is incorporated in luminescence 
markers and is being considered for use in nuclear batteries. 

Both of the standard samples may be ordered under the general 
licensing provisions of the Atomic Energy Act of 1954. The iron 55 
standard sample No. 4929 and the promethium 147 standard sample 
No. 4940 may be purchased for $36.00 and $40.00, respectively. Both 
samples are in solution form (approximately 3 ml), and have nominal 
activities of 5 X 104 dps/g. The iron 55 solution contains lug/ml 
FeCl; as carrier, while the promethium 147 solution is carrier free. 
All orders for the radioactivity standards should be addressed to (Miss) 
Elizabeth M. Zandonini, Radioactivity Standard Samples, Radio- 
activity Section, National Bureau of Standards, Washington 25, D. C. 
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THEORY OF MACHINES, by Joseph E. Shigley. 
Part I, KINEMATIC ANALYSIS OF MECHA- 
NISMS, 631 pages, diagrams, illustrations, 
6 X 9h in. New York, McGraw-Hill Book 
Co., Inc., 1961. Price, $11.50. Parr II, 
DyNAMIC ANALYSIS OF MACHINES, 645 
pages, diagrams, illustrations, 6 X 9} in. 
New York, McGraw-Hill Book Co., [nc., 
1961. Price, $7.50. 


This very complete and well-written book 
treats the standard as well as the more ad- 
vanced topics of machine analysis with great 
clarity. The use of vector notation is intro- 
duced in the first chapter to lay a solid 
foundation for the exposition of kinematics, 
and of the geometry of motion, which com- 
prise over one third of Part 1. This is fol- 
lowed by chapters on gearing, gear trains, 
linkages and cams which, in keeping with the 
purpose of the book, are discussed from the 
point of view of analysis rather than design. 
Many references are given to more detailed 
design applications. The advantages of 
vector methods are clearly illustrated in the 
chapter on space mechanisms. A brief dis- 
cussion of analog computing mechanisms 
concludes Part 1. 

Part 2 applies and adapts the more ad- 
vanced concepts and techniques of dynamic 
analysis to a representative selection of prac- 
tical problems. Vector methods are con- 
sistently used to solve dynamic problems in 
plane and space mechanisms. The discussion 
of the phase-plane delta method for analyzing 
transient vibration phenomena is clear and 
comprehensive. Considerable space is de- 
voted to the dynamics of rotating machinery, 
including problems of balancing, and to cam 
mechanisms composed of non-rigid com- 
ponents. A brief introduction to the dy- 
namics of feedback control systems concludes 
the book. 

Comparable works have seldom satisfied, 
at the same time, the rigorous standards of 
the theoretician and the demands of the 
practicing engineer, for techniques applicable 
to a great variety of technical problems. 
Professor Shigley has succeeded in presenting 
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the analytical foundations of machine design 
in a manner simple enough to be of value to 
the student, yet so comprehensive and general 
that it may serve as a reference to the design 
engineer. 

KLaus CAPPEL 
The Franklin Institute Laboratories 


LECTURES ON LINEAR ALGEBRA, by I. M. 
Gel’fand. 185 pages, 6X9 in. New 
York, Interscience Publishers, Inc., 1961. 
Price, $6.00. 


The last few years have seen the publica- 
tion of a good many books bearing the same 
or a very similar title. These books vary 
widely in their approach and their style. 
Many are written in such an abstract and 
“uncompromising”? style that only pro- 
fessional mathematicians can read them with 
profit. 

The merit and charm of this book lie in 
the rare blend of mathematical elegance, 
pedagogy and compactness. The author 
prepares the reader to the next development 
and motivates the extensions and definitions 
before stating them. He frequently appeals 
to intuitions and analogies, leading naturally 
to generalized concepts. No mathematician 
could question the rigor of the theorems, 
proofs and definitions as they are exposed, 
yet the reader has a feeling of participation 
in the build-up of the theory. Carefully 
chosen and very diverse illustrative examples 
contribute to improve the understanding and 
dispel possible misgivings. A surprising 
amount of material is covered in less than 
two hundred pages. Let us mention in 
particular: real and complex n-dimensional 
spaces; Euclidean spaces; reduction of quad- 
ratic forms (a very thorough treatment) ; 
linear, hermitean and unitary transforma- 
tions; the eigenvalue problem (extremal 
properties) and an introduction to the algebra 
of tensors. 

The book ends rather abruptly and the 
chapter on tensors seems unfinished. Two 
important appendices have been deleted (at 
the author’s request) from the English edition, 
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and the publishers (in this otherwise carefully 
printed and bound book) have omitted an 
index of topics. Forgetting these 
defects, the book appears to be a very suitable 
introduction to this important branch of 
mathematics and mathematical physics. 
HENRI AMAR 
Temple University 
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FLow oF FLurps THRouGH Porous MATE- 
Royal E. Collins. 
diagrams, illustrations, 6 X 9 in. 
York, Reinhold Publishing Corp., 
Price, $12.50. 


263 pages, 
New 
1961. 


RIALS, by 


In the preface of this text the author 
indicates the importance of flow through 
porous materials in many phases of petroleum 


ground-water 
This 


engineering, soil mechanics, 
hydrology and sanitary engineering. 
book is intended as a text for the petroleum 
reservoir engineer specifically. The format 
consists of ten chapters dealing with struc- 
ture, statics, physical and mathematical 
theory, steady and transient laminar flow 
of homogeneous fluids, flow 
of immiscible fluids, moving boundary prob- 


lems, theory of models and flow with phase 


simultaneous 


change. 

The treatment of each chapter is a highly 
theoretical, mathematical development be- 
yond the capabilities of the average under- 
Every re- 

technical 


graduate scientist or engineer. 
viewer has favorite portions of 
publications and in this case it was the 
chapters concerned with moving boundary 
problems and flow with change of phase. 
The former should be of interest to petroleum 
producing engineers, deep-well water investi- 
gators and others dealing with similar prob- 
lems of calculating complex drawdown pat- 
terns in multi withdrawal fields. Less obvious 
is the extension of the theoretic concepts to 
the complex drying or extraction process con- 
trolled by unsteady state parameters. Ad- 
mittedly, the latter chapter is treated rather 
briefly, but the discussions of phase distribu- 
tion and saturation will be of interest not 
only to the obvious audience, but to those 
researchers in chromatography and solvent 
extraction. 

The many derivations are smoothly pre- 
sented with adequate explanatory wording to 
permit the practicing researcher or student 
to follow without undue cross-referencing or 
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assumptions. This is an excellent treatment 
of a subject that has not been too well 
elaborated by the writers. Al- 
though it describes basic concepts applicable 
to many fields, the advanced and theoretical 
nature of the text will limit its circulation to 
research oriented readers. 
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Ropert A. BAKER 
The Franklin Institute Laboratories 


Tue THeory oF TRANSITION-METAL Ions, 
by J. S. Griffith. 
7X10 in. New York, Cambridge Uni- 
versity Press, 1961. Price, $17.50. 


455 pages, diagrams, 


The subject matter of this excellent work 
is somewhat less general than the title would 
indicate, since it is restricted to what is called 
ligand-field theory. Metallic compounds and 
alloys of transition elements, and, of course, 
the pure elements themselves, are excluded. 
The theory applies to situations in which the 
states which accommodate the d- or f-orbitals 
of the are, 
to a good approximation, localized on the 


isolated transition element ion 
central transition element ion and its nearest 


neighbors or ligands. Cooperative phe- 
nomena involving several transition element 
ions are also excluded from the treatment. 

After a brief introductory chapter, the 
author begins with the quantum mechanical 
description of the angular momentum of the 
free atom. The bra and ket notation of 
Dirac is used here and throughout the book. 
There is a short chapter on electromagnetic 
radiation and its interaction with an isolated 
Chapters 4 and 5 with the 
structure of free atoms and ions and with 
magnetic effects, respectively. Spin is in 
troduced first phenomenologically and then 
via Dirac’s wave equation. Chapter 6 con- 
cludes the survey of background material 
with a treatment of group theory and matrix 
representations of groups. An effort is made 
to collect useful theorems and formulae for 
the benefit of those using the work as a refer- 
ence only. These collections are unfortu- 
nately not easily found without paging through 
the text. 

With Chapter 7 the concept of complex ion 
is introduced. ‘The possible coupling schemes 
are discussed. These depend on the order 
in which the electrostatic interaction between 
the partially filled shells and the rest of the 
ion (U), the mutual electrostatic interaction 
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of the electrons in the ion-filled shells (3C,), 
and the spin-orbit interaction energy (J5C,) 
are diagonalized. If the order is K,, 3,, VU, 
the treatment is equivalent to crystal field 
theory. ‘This order is termed the rare earth 
coupling scheme for obvious reasons. In 
the weak field coupling scheme the order is 
UL, and in the strong field coupling 
scheme it is 5C,, The first two schemes 
are treated in Chapter 8 and the third in 
Chapter 9. In Chapter 10 the paramagnetic 
susceptibility of different electronic con- 
figurations are calculated, and it is shown 
how the experimental determination of sus- 
ceptibility often leads to an unambiguous 
spin assignment to the electronic configura- 
tion of the complex ion. Optical spectra 
and thermodynamic properties of some com- 
plex ions are calculated and compared with 
experimental results in Chapter 11. Chapter 
12 is devoted to the calculation of the electron 
and nuclear spin resonance spectra of the com- 
plex ions. Appendices are provided for special 
topics not treated in the text. 

The book is well written and referenced. 
The explanations are clear and the style, 
after the first shock of meeting the first 
person singular, stimulating. Exercises are 
provided after most of the sections. This 
is a valuable textbook and useful reference. 

F. J. DONAHOE 
The Franklin Institute Laboratories 


A HOLE IN THE BoTTOM OF THE SEA, by 


Willard Bascom. 352 pages, diagrams, 
illustrations, plates, 53 X9} in. New 
York, Doubleday & Co., Inc., 1961. Price, 


$4.95. 


The statement that we know more about 
the distant stars than we do about the earth 
on which we live is an obvious truth. Analy- 
tical boring machines have yielded a fairly 
comprehensive picture of the inside of a star. 
Not to be outdone, the earth scientists are 
trying to unlock the secrets of the earth to 
provide more knowledge of our planet. They 
have conjured up a fantastic scheme proposed 
by as unlikely a society as can be imagined 
and the combination of the two has resulted 
in this country being involved in the drilling 
of a hole in the bottom of the sea. 

For a long time, seismologists have tried 
to determine the structure of the earth. 
Shrewd guesses and skillful interpretations 


Book REVIEWS 


331 


of seismic waves have yielded some knowledge 
of the earth’s structure. In the interpreta- 
tions, scientists discovered a discontinuity 
which comes within 6 miles of the sea’s sur- 
face. Because the earth is so much thicker 
on land, the obvious approach to pierce the 
surface layer is to drill the hole under the 
sea. Thus results the attempt to bore 
through the upper crust of the earth to un- 
cover the composition of the subcrustal 
layers. There is a bonus to be taped in this 
process. When the drill brings up cores of 
the ocean sediments it may prove to be a key 
to millions of years of the past life history 
of the earth. 

A Hole in the Bottom of the Sea is the story 
of the origin and the implementation of the 
drilling of the Mohole. It is the story of the 
scientists’ attempt to breach the crustal 
layers of the earth to learn more about the 
structure. In a magnificently written book, 
Willard Bascom has chronicled the story of 
how the idea originated and what is being 
done and what scientists hope to do. The 
author tells his story in a fascinating manner 
and skillfully weaves a tremendous volume 
of information in the course of telling the 
story. Even the technical details in the art 
of drilling holes both on land and on sea prove 
excellent exposition. The use of line draw- 
ings proves a valuable asset in visualizing the 
rigs and their uses. 

Despite the fact that the author is the 
director of the Mohole Project, he does not 
paint too rosy a picture concerning the suc- 
cess of the project. He indicates the tremendous 
technical difficulties which will be encoun- 
tered and possible solutions. Fortunately 
he is a realist and recognizes fully the mag- 
nitude of the task. 

Several minor errors have crept into this 
work when the author gets involved in as- 
tronomy. He speaks of the shift of the red 
calcium band. The calcium lines are in the 
blue. He speaks about the farthest stars 
moving fastest away. The stars do not move 
in the sense of all moving away from the earth 
—it is the star systems or galaxies to which 
the author refers. In his table of planets 
and moon the number of Jovian moons is 
given as 11—there are 12. 

These are minor items in one of the most 
fascinating books it has been this reviewer's 
pleasure to read. For those who are inquisi- 
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tive about one of the more startling adven- 

tures of science they can do no better than 

to read the exciting A Hole in the Bottom of 
the Sea. 

I. M. Levitt 

The Fels Planetarium 


THe DesiGN oF CHEMICAL REAC- 
rors. VoL. 3, MATHEMATICS IN SCIENCI 
AND ENGINEERING, by Rutherford Aris. 
191 pages, diagrams, 6 KX 9in. New York, 
Academic Press Inc., 1961. Price, $7.00. 


Advanced mathematical concepts coupled 
up with computer solutions have become in- 
creasingly important in the solution of com- 
plex engineering problems. Dynamic _pro- 
gramming has proven especially powerful 
and versatile in solution of the engineering 
rhe 
practicing design engineer and mathematician 
have had to cooperate to a greater extent 
than ever before, but often without complete 
comprehension of the other’s language. 
Elucidation of this common language has 
been effectively accomplished for both by 


and economic aspects of reactor design. 


this treatise. 

The author's personality comes through 
in a manner not common in scientific dis- 
cussions of depth. Chapters are prefaced 
by witty rhymes, an excellent style with 
choice the subject 
matter, and there are liberal cautions against 
building a mathematical house of cards based 
The content 


vocabulary describes 


on unnecessary complications. 
is very mathematical as befits the subject. 
The reader should be grounded in such basic 
advanced mathematical concepts as proba- 
bility distributions, transforms, etc., to follow 
the derivations. An understanding of chem- 
ical engineering kinetics and thermodynamics 
will also be useful. 

Reactor designs have previously been de- 
veloped through trial and error experience. 
As understanding has evolved, basic concepts 
have been developed in mathematical format. 
The modern reactor design problem involves 
many parameters and manpower and time 
considerations preclude the classical engineer- 
ing approach. Reduction of the overall 
problem to manageable portions for solution 
by digital computers via dynamic program- 


ming has proven ideal. Since the design 
premises may be stated with fairly high 
accuracy, engineering design by dynamic 
methods is sounder than solution in other 
fields such as economics. 

There are three components of optimal 
design. The specification of: (1) process 
stream state, (2) operating variables and 
transformations they effect, and (3) objective 
function to be optimized. In this book 


functional equations suitable for digital com- 
puter solution are evolved for a variety of 


reactor types. Consideration is given to 
physical and chemical factors with the num- 
ber of independent reactions being indicative 
of the problem complexity. Cost estimation 
procedures are evolved based on (1) capital 
and maintenance costs (2) stream costs (3) 
heat costs and (4) control costs. Functions 
are developed for unit profit per unit time. 
Optimal equations are derived for continuous 
multibed adiabatic 
Temperature 


flow stirred reactors, 
reactors and tubular reactors. 
and holding time, feed stream bypassing, 
interchange and cold shot cooling, etc. are 
considered where applicable. Stochastic ap- 
plications are given and two in particular 
show the treatment. These 
are optimal replication of processes subject 
to failure, and optimal catalyst replacement 
policies. In the latter case a management 
policy is obtained through maximization of 
catalysts replacement cost. A 


need for such 


profit vs. 
brief treatment of optimizing operation of 
existing reactors completes the technical 
content divided into nine chapters. 
References are given at the end of the book 
and it is vividly clear that only a handful 
of authors all of current vintage are repre- 
sented. This is a new, important subject. 
Chemists, chemical engineers, designers, oper- 
ations analysts and computer specialists will 
benefit by this work. The publisher’s cover 
includes undergraduates in the suggested 
audience. Unless the undergraduate educa- 
tional processes are suddenly accelerated with 
computer speed this is likely to be an over- 
Scientists and math- 
valuable 


enthusiastic statement. 
ematicians will find this a 
addition to their libraries. 


very 


R. A. BAKER 
The Franklin Institute Laboratories 
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REFERENCE ELECTRODES: THEORY AND PRAC- 


TICE, edited by David J. G. Ives and 
George J. Janz. 651 pages, diagrams, 
6X9 in. New York, Academic Press, 
1961. Price, $20.00. 


This new graduate-level text is principally 
concerned with aspects of 
electrochemical theory. 
are examined from both the mechanistic and 
separation 


experimental 
Reference electrodes 


thermodynamic viewpoints. 
is made between theory and experiment; 
each chapter contains both. After an intro- 
ductory chapter, eight chapters deal with 
specific kinds of electrodes—for example, 
hydrogen, calomel, silver halide, quinhydrone 
—and the three final chapters cover general 
topics such as reference electrodes in non- 
aqueous solutions, biological uses of micro- 
electrodes, and electrodes in fused salt sys- 
tems. Not the least attractive feature of 
the book is the extensive list of references ap- 
pended to each chapter. The volume, in 
addition to its usefulness in filling a gap in 
graduate texts, will be of interest to research 
workers. 


Puysics, VoL, 
Brittin, B. W. 
531 pages, 


LECTURES IN THEORETICAL 
III, edited by Wesley E. 
Downs and Joanne Downs. 


diagrams, 6 X 9 in. New York, Inter- 
science Publishers, Inc., 1961. Price, 
$11.00. 


Volume III of this series contains lectures 
delivered at the Summer Institute for Theo- 
retical Physics at the University of Colorado 
in 1960. The level of presentation of the 
eleven lectures lies between that of the text- 
book and that of the research journal. The 
wide variety of subjects and the eminence of 
the contributors are evident from this listing 
of the papers: “Causality and Dispersion 
Relations,” by Aage Bohr; ‘‘Selected Topics 
in Theoretical Physics,” by Victor Weisskopf ; 
“Statistical Mechanics of Irreversibility,” by 
Robert Zwanzig; ‘The Scattering of Elec- 
trons by Atoms,’”’ by B. L. Moiseiwitsch; 
“Some Applications of the Generating Func- 


BOOK NOTES 


tional of the Molecular Distribution Func- 
tions,” by Melville Green; ‘Nonequilibrium 
Statistical Mechanics,” by Elliott Montroll; 
“Quantum Theory of Collision Processes,” 
by R. Haag; “Canonical Commutation Rela- 
tions in Field Theory and Functional Inte- 
gration,” by R. Haag; ‘‘General Theory of 
Nonequilibrium Phenomena,” by Radu Ba- 
lescu; ‘“‘Martin-Schwinger Methods in the 
Many-Body Problem,” by Neil Ashby; and 
“Green's Functions and the Quantum Theory 
of Fields,”” by Kurt Symanzik. 


INTRODUCTION TO ELEMENTARY PARTICLE 
Puysics, by R. E. Marshak and E. C. G. 
Sudarshan. Interscience Tracts on Phys- 
ics and Astronomy #11, 217 pages, tables, 


5} X 8 in. New York, Interscience Pub- 
lishers, 1961. Price: $2.50 (paper); $4.50 
(cloth). 


This tract concentrates on those aspects 
of elementary particle physics which do not 
depend on detailed dynamical calculations. 
The physical meaning of the various concepts 
is emphasized and the experimental material 
is introduced chiefly for illustrative purposes. 
The book is intended as a brief introduction 
for the non-expert in theoretical high energy 
physics. 


PHYSICAL CHEMISTRY, by Gordon M. Barrow. 
684 pages, diagrams, illustrations, 6} X 9} 
in. New York, McGraw-Hill Book Co., 
Inc., 1961. Price, $8.95. 


An introductory, undergraduate text, this 
work deals with all the subjects customarily 
dealt with, such as thermodynamics, states 
of matter and kinetics. The author unifies 
the various subjects by interpretations based 
on molecular structure and behavior, for 
he believes that physical chemistry is not a 
branch of chemistry, but rather a method 
which is applicable to all branches of chemis- 
try. The text is roughly divided into two 
parts; the first deals with situations involving 
a single or a few molecules at a time, while 
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the second involves the interplay of many 
molecules. 


HANDBOOK OF THERMOPHYSICAL PROPERTIES 
OF SOLID MATERIALS, compiled by Alex- 
ander Goldsmith, Thomas E. Waterman 

and Harry J. Hirschhorn. Five volumes, 

revised edition, 4300 -pages, diagrams, 
charts, 8$ X11 in. New York, The 

Macmillan Co., 1961. Price, $90.00 (five 

volumes). 


Prepared by Armour Research Foundation 
under contract to the U. S. Air Force, this 
five-volume reference work provides a much 
needed source of reliable thermophysical 
property data on high temperature materials. 
Twelve properties are included: 
melting point, latent heat of fusion, latent 
heat of vaporization, latent heat of sublima- 
tion, specific heat, thermal conductivity, 
thermal diffusivity, emissivity and_ reflec- 
tivity, linear thermal expansion, vapor pres- 
The first four 


density, 


sure, and electric resistivity. 
volumes deal, respectively, with: [, elements 
(melting points above 1000° F.); II, alloys; 
IIT, ceramics; and [V, cermets, intermetallics, 
polymerics, and Volume \ 
contains the list of references from which the 


composites. 


data were compiled and an author index. 
The data are presented clearly, both in 
tabular and graphical form. Workers in 
missiles, aircraft and the nuclear sciences will 
find this handbook invaluable and easy to 
use; the sturdy binding will withstand hard 


usage. 


AN INDEX OF PUBLISHED INFRA-RED SpPEC- 
rRA, Vous. I ANb II, edited by Ministry of 
Technical Information and Li- 
brary Services. 805 pages, 8} X 13} in. 
London, Her Majesty’s Stationery Office, 
1960. Price, $18.40. American 
British Information Services. 


Aviation 


agents, 


This ‘Beilstein’ type of index enables a 
given spectrum to be traced in the literature. 
The first two volumes include most of the 
spectra published up to 1957, and subsequent 
volumes will follow. All the infra-red spectra 
are listed under a formula index; the state 
in which the sample was studied is also given, 
as are the wavelength range covered and an 
indication of the optical system employed. 
Nearly 10,000 references are covered. Prac- 
ticing spectroscopists will find this a valuable 
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source of information, and a real time- 


Saver. 


BLack GoLp, by Arthur Beeby-Thompson. 
533 pages, plates, 64 X 9} in. New York, 
Doubleday & Co., Inc., 1961. Price, $5.95. 


This is the autobiography of the Grand 
Old Man of oil, covering the exciting early 
days of this industry, beginning in 1898, up 
through the latest techniques. The reader 
is taken to Africa, Baku, the Caribbean, 
Mesopotamia, Peru, Egypt, America and 
Arabia, as the author recounts his adventures 
in pioneering in the oil industry. 


MANAGEMENT AND INDUSTRIAL 
APPLICATIONS OF LINEAR PROGRAMMING, 
Volume 1, by A. Charnes and W. W. 
Cooper. 467 pages, diagrams, 7 X 104 in. 
New York, John Wiley & Sons, Inc., 1961. 
Price, $12.50. 

The contents of this new work are based 
on research conducted by the authors for 
almost ten years. It is directed mainly to 
persons who are interested in managerial 
applications of linear programming, and is 
intended to serve as a textbook as well as a 
reference. 

Some of the general topics covered are 
the Simplex Method and Some Geometry 
in Linear Spaces, Data Interpretations by 
Direct and Dual Horizons and 
Surrogates in a Dynamic Model for Produc- 
tion Scheduling, and Degeneracy, Perturba- 
tion, and Complete Regularization. 


Relations, 


QuANTUM THEORY: ELEMENTS, VOLUME 1, 


edited by D. R. Bates. 435 pages, dia- 
grams, 6 X9 in. New York, Academic 
Press Inc., 1961. Price, $10.00. 


he first in a series of three volumes, this 
book is intended as an advanced text and 
reference on the fundamentals and applica- 
tions of quantum theory. The reader should 
have a knowledge of classical mechanics, 
electromagnetic theory and atomic physics, 
and be familiar with the ordinary techniques 
of mathematical analysis. 

In this first volume of Quantum Theory, 
non-relativistic wave mechanics and matrix 
mechanics are introduced and an extensive 
survey of the exactly soluble problems of 
the point and of the continuous spectrum is 
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presented. Developed in considerable detail 
are the approximate methods which are 
available for treating other stationary and 
time-dependent problems. Also treated fully 
is scattering theory, meeting the needs of 
both ionic and nuclear physicists. 


OPERATIONAL ELeEctricITty, by Charles I. 
Hubert. 510 pages, diagrams, illustra- 
tions, 6 X 9 in. New York, John Wiley & 
Sons, Inc., 1961. Price, $8.50. 


Intended as a text for students enrolled 
in a two-semester nonelectrical engineering 
course, this book deals with theory, character- 
istics, applications, and mode of operation of 
circuits and machines. 

The author, who has had 18 years teaching 
experience in electrical engineering and is an 
associate professor of that course at the 
U. S. Merchant Marine Academy, has in- 
cluded in his latest book chapters on Circuit 
Elements, Balanced Polyphase System, Di- 
rect-Current Generators, Alternating-Current 
and Direct-Current Motors, Generator Volt- 


age Regulators, Industrial Electronics and 
Instrumentation. 
PROGRESS IN VERY HIGH PRESSURE RE- 


SEARCH, edited by F. P. Bundy, W. R. 
Hibbard, Jr., and H. M. Strong. 314 
pages, diagrams, illustrations, 8 X 11} in. 
New York, John Wiley & Sons, Inc., 1961. 
Price, $12.00. 


This book contains the proceedings of an 
International Conference held on June 13-14, 
1960, at Bolton Landing, Lake George, New 
York. The conference was sponsored by the 
Materials Central, Wright Air Development 
Division of the U. S. Air Force and the 
General Electric Research Laboratory. 

Since, in recent years, there has been a 
considerable increase in the use of very high 
pressure as a research method, it was felt 
that much progress could be achieved by 
bringing together some of the most active 
workers in the field for exchange of data, 
experience and philosophy. editors 
hope that this book will help to stimulate and 
guide new work, and to encourage the con- 
tinuation of current investigations. 

The central topics of the conference were 
equipment and technique, structures of 
materials formed under high pressure, and 
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behavior of matter under high pressure. 
Of particular interest in this book is a review 
of high pressure work in the USSR, presented 
by Professor L. F. Vereshchagin of Moscow. 


TEMPERATURE MEASUREMENT IN ENGINEER- 
ING, Vor. II, by H. Dean Baker, N. H. 
Baker and E. A. Ryder. 465 pages, 
diagrams, 6 X9 in. New York, John 
Wiley & Sons, Inc., 1961. Price, $13.50. 


Volume II of this work completes the task 
begun in Volume I, published in 1953. Vol- 
ume I dealt with thermocouples; Volume 
II is devoted to resistance thermometry 
and radiation pyrometry, as techniques in 
the measurement of temperature. Both are 
treated descriptively and analytically, as 
general techniques, in the first part of the 
book. The second, and larger, part contains 
a selected assortment of problems chosen 
from among the almost infinite variety 
available. In the applications portion of the 
book, chapters are devoted to such fields 
as surface temperature, moving bodies, 
liquids, gases, nuclear-reactor temperatures, 
and flame temperatures. In these, concrete 
examples of successful techniques are given; 
few worked-out examples are included. 


THE CoL_LtecteD Works OF IRVING LANG- 
MUIR, compiled under the general editor- 
ship of Guy Suits. Twelve volumes, 
varying from approximately 250 to 500 
pages, diagrams, 7 X 9} in. Oxford, Per- 
gamon Press, 1960 & 1961. Price, $15.00 
(£ 5), each volume. 


This twelve-volume edition of the collected 
works of the first American industrial scien- 
tist to receive a Nobel Prize is being published 
under the general editorship of Guy Suits, 
one of Langmuir’s associates in the General 
Electric Company’s Research Laboratory. 
Langmuir’s career was so varied and his 
scientific output so prolific that the editors 
have divided his contributions by subject 
matter, for inclusion in separate volumes. 
The first six, copyrighted in 1960 and 1961, 
are now available. They deal, respectively, 
with: 1, Low-Pressure Phenomena; 2, Heat 
Transfer—Incandescent Tungsten; 3, Therm- 
ionic Phenomena; 4, Electrical Discharge; 5, 
Plasma and Oscillations; and 6, Structure of 
Each volume contains, in addition 
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to Langmuir’s own papers, an article analyz- 
ing his contributions to the particular field 
covered in that volume. 

The editors and the publisher are to be 
congratulated for the care with which this 
important scientist’s widely scattered con- 
tributions have been brought together for 
systematic presentation. 


THE CELL AND THE ORGANISM, edited by 
J. A. Ramsay and V. B. Wigglesworth. 
350 pages, diagrams, plates, 6 X 9} in. 
New York, Cambridge University Press, 
1961. Price, $9.50. 

This is a collection of essays on the applica- 
tion of experimental methods to different 
types of investigation in the zoological field. 
The subject matter covers a wide range, from 
the study of individual cells to the mechanics 
of motion in vertebrates. The contributions, 
written in honor of Sir James Gray, one of 
zoology’s most eminent men, are original and, 
in some cases, speculative; they make a 
valuable addition to the literature on zoology 
and biology. 

ENGINEERING MANAGEMENT AND ADMINIS- 

TRATION, by Val Cronstedt. 335 pages, 

diagrams, 6 X 8 in. New York, McGraw- 

Hill Book Co., Inc., 1961. Price, $8.50. 


Administrators who find themselves in 
charge of an engineering operation will find 
much of interest in this new book. Covering 
all phases of engineering administration 
except research, the book contains many 
helpful suggestions (on budget, proposals, 
test operations, patent law, records, etc.), 
designed to free the technical man from a 


military type of discipline. 


[HE CILIATED Protozoa, by John O. Corliss. 
Volume 7 of the International Series of 
Monographs on Pure and Applied Biology, 

Division, 310 pages, plates, 
Oxford, Pergamon 

Price, $12.00. 


This monograph is written at the graduate 
level, but it should also be useful to teachers, 
post-graduate students and research workers 
in protozoology, parasitology and cell biology. 
The author proposes a new scheme of classifi- 
cation of the ciliated protozoa, based on new 
facts and revised interpretation of older 


Zoology 
53 X 82 in. 
Ltd., 1961. 


Press 
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data. Since this is the first work devoted 
solely to ciliated protozoa since 1935, such 
a reclassification seems called for, so that 
new families, genera and species may be 
fitted into the classification. A systematic 
over 1400 scientific 
over 1600 


names of 
literature 


index lists 


ciliate groups, and 


references are cited. 


MECHANICAL ELE- 
509 pages, 
New 
1961. 


OpTIMUM DESIGN OF 
MENTS, by Ray C. Johnson. 
diagrams, illustrations, 6 X 9 in. 
York, John Wiley & Sons, Inc., 
Price, $11.50. 


The author presents a method of mechan- 
ical design, a result of his many years of 
industrial and academic experience, which 
is explicit in nature and can be applied to 
innumerable original design problems. His 
method minimizes the “‘cut-and-try’’ prac- 
tice, and takes into account the inherently 
unavoidable limitations which confront prac- 
ticing mechanical design engineers. 

rhe thirteen chapters include: Effects of 
Manufacturing Errors on Product Perform- 
ance, Optimum Choice for Method of Analy- 
sis, Optimum Design of Torsion Shafts, and 
Optimum Design of Gears. 


LINEAR SYSTEMS ANALysIs, by Paul E. 
Pfeiffer. 538 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1961. Price, $12.50. 


Intended to provide a theoretical core of 
fundamentals needed to understand a variety 
of systems and devices, this text stresses 
tools and habits of analytical 
thought. The author this through 
stressing: (1) physical analogies; (2) general 
self- 


analytical 
does 


rather than specific cases; and (3) 
contained mathematical treatment. As a 
text, the first seven chapters can be used as 
a one-semester introductory course for 
juniors; the entire book can be used on the 


senior level. 


ELECTROMECHANICAL SYSTEM THEORY, by 
Herman E. Koenig and William A. Black- 
well. 497 pages, diagrams, 6 X9 in. 
New York, McGraw-Hill Book Co., Inc., 
1961. Price, $14.50. 


The authors have prepared this text to 
present a new engineering systems discipline 
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which encompasses the wide variety of com- 
ponents found in the rapidly expanding con- 
trol-systems field. In contrast to the notions 
of analogies used for many years, the applica- 
tion of linear graph theory is stressed as the 
core of the authors’ proposed discipline. In 
this text, the subjects of rotating machinery, 
polyphase and servomechanisms 
have been integrated into the discipline, as 
well as those aspects of mechanics, hydraulics 
and electronic circuitry which are applicable 
Emphasis is placed on 


systems 


in control systems. 
precise formulation procedures inherent in 
linear graph theory. Prerequisites for study 
of the text are matrix algebra and Laplace 
transform techniques. 


ADVANCES IN ELECTRONICS AND ELECTRON 
Puysics, Vou. XIV, edited by L. Marton 
and Claire Marton. 329 pages, diagrams, 
illustrations, 6 X 9 in. New York, Aca- 
demic Press Inc., 1961. Price, $11.00. 
The fourteenth volume of this well-known 

series contains the following six contributions: 

“The Electron as a Chemical Entity,” by 

C. G. B. Garrett; “Problems of Photocon- 


ductivity,”’ by P. Gérlich; ‘Strong-Focusing 


Septier; ‘Hydrogen 
Goldberg and 


Albert 
Seymour 


Lenses,” by 
Thyratrons,” by 
Jerome Rothstein; “Cerenkov Radiation at 
Microwave Frequencies,”’ by Herbert Lashin- 
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sky; and “High-Power Axial-Beam Tubes,’ 
by T. Moreno. 


HANDBOOK OF SUPERSONIC AERODYNAMICS. 
SEcTION 20. WIND TUNNEL INSTRUMEN- 
TATION AND OPERATION, by R. J. Volluz. 
511 pages, diagrams, illustrations, 74 X 10} 
in. Washington, Government Printing 
Office, 1961. Price, $3.00. 

Section 20 of this Handbook compiled at 
Johns Hopkins is part of Volume 6, which 
covers design, construction and operation 
of the test facilities (free jets, shock tubes and 
wind tunnels) from which are obtained aero- 
dynamic data. Section 20 deals with wind 
tunnel instrumentation and operation (de- 
sign will be covered in Section 19). Three 
chapters cover in detail: various methods 
for measuring pressure, force, and tempera- 
ture; the effects on such measurements 
produced by yaw, density, load, etc.; and 
errors associated with such measurements. 
Other chapters deal extensively with flow 
(Schlieren systems, interferometry, boundary- 
layer transitions, etc.), turbulence, tunnel 
blockage, model design, typical force measur- 
ing systems (reduction and plotting tech- 
niques), and pressure measuring systems 
(pressure transducers, manometers, etc.). 
This Section of the Handbook maintains the 
high standards set in the previously published 
sections. 


Set 
i 
2 
: 
+ 
| 


PUBLICATIONS RECEIVED 


FreLp GroLoGy, by Frederic H. Lahee. Sixth edition, 926 pages, diagrams, illustrations, 
5 X 74in. New York, McGraw-Hill Book Co., Inc., 1961. Price, $10.75. 


REACTOR HANDBOOK, VOL. IT, FUEL REPROCESSING, edited by S. M. Stoller and R. B. Richards. 
Second edition, revised and enlarged, 665 pages, diagrams, illustrations, 7 XK 10 in. New 
York, Interscience Publishers, Inc., 1961. Price, $21.40. 


THE Story OF ATOMIC THEORY AND ATOMIC ENERGY (Formerly titled THe Atom Srory), 
by J. G. Feinberg. Unabridged republication of the first edition, with a new Foreword 
and an added chapter, 264 pages, diagrams, 5} K 8 in. New York, Dover Publications, 
1961. Price, $1.45 (paper). 


Space & Time, by Emile Borel. Unabridged and unaltered republication of the first English 
translation, with a new Foreword by Banesh Hoffman, 234 pages, 5} K 8 in. New York, 
Dover Publications, 1961. Price, $1.45 (paper). 


PIONEERS OF SCIENCE, by Sir Oliver Lodge. Unabridged and unaltered republication of the 
last corrected edition published in 1926. 404 pages, illustrations, 5} & 8in. New York, 
Dover Publications, 1961. Price, $1.50 (paper). 


WATER HAMMER IN HYDRAULICS AND WAVE SURGES IN ELecrricity, by Louis Bergeron. 
Translated from the 1949 French edition under the sponsorship of the Am. Soc. of 
Mech. Engineers, 293 pages, diagrams, inserts, 83 X 11 in. New York, John Wiley & 
Sons, Inc., 1961. Price, $15.00. 


FUNDAMENTALS OF PHysICAL GEOGRAPHY, by Glenn T. Trewartha, Arthur H. Robinson and 
Edwin H. Hammond. 389 pages, diagrams, illustrations, maps, 7 X 9} in. New York, 
McGraw-Hill Book Co., Inc., 1961. Price, $6.95. 


SoME Aspects OF LIFE IN FRESH WATER, by Edward J. Popham. Revised edition, 127 pages, 
diagrams, 4} X 7} in. Cambridge, Harvard University Press, 1961. Price, $2.25. 


NEw SuHortcuts To tv SERVICING, by Leonard C. Lane. Two volumes, 160 pages each, 
diagrams, 54 X 83 in. New York, Gernsback Library, Inc., 1961. Price: $3.20 each 
volume or $5.90 for both volumes (paper) ; $9.90 for both volumes (cloth). 


A To Z 1N Aupto, by G. A. Briggs. 224 pages, diagrams, 54 84 in. New York, Gernsback 
Library Inc., 1961. Price: $3.20 (paper): $5.00 (cloth). 


SCIENCE AND THE NEW Nations, edited by Ruth Gruber. 314 pages, plates, 6 X 9} in. 
New York, Basic Books, Inc., 1961. Price, $6.50. 


QUANTITATIVE ORGANIC MICROANALYsIS, by Al Steyermark. Second edition, revised and 
enlarged, 665 pages, diagrams, illustrations, 6 X 9 in. New York, Academic Press Inc., 
1961. Price, $16.50. 
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Sapphire Whiskers for Silver.— 
Space scientists have increased the 
strength of pure silver five times by 
blending into it a small percentage of 
tiny sapphire ‘‘whiskers.”” At the 
same time they settled a long-standing 
question of whether the ‘‘whisker”’ 
technique can be used to produce 
ultra-strength metal materials. 

Dr. Willard H. Sutton, of General 
Electric’s Missile and Space Vehicle 
Department in Philadelphia, said 
similar techniques applied to high- 
strength metals can be expected to pro- 
duce materials several times stronger 
at white heat than any high-strength 
alloy. 

The purpose of the work, being 
done at the Department’s Space 
Sciences Laboratory under a contract 
from the Bureau of Naval Weapons, 
is to help overcome one of the prob- 
lems in space vehicle and advanced air- 
craft design—the heat barrier. There 
is a lack of structural materials that 
will maintain strength at the high 
temperatures created by reentry, hy- 
personic air speeds, and powerful jet 
and rocket engines. 

The theory behind ‘‘whisker’’ rein- 
forcement is: 

Laboratory-grown whiskers, which 
consist of a single crystal of alumina 
(or sapphire), have shown extremely 
high strength—some can withstand 
tensions up to several million pounds 
per square inch. However, they have 


neither the rigidity nor the bulk to 
be machined or used as structural 
pieces or engine parts. 

Currently known high-strength met- 
als, including steels, can be machined 
into precision parts and are rigid, but 
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at elevated temperatures (2000 deg. 
F. or more) they lose nearly all of 
their strength. 

By embedding the whiskers inside 
the metal properly, the whiskers will 
carry nearly all of the load, thus 
permitting a metal to perform at 
temperatures far beyond its normal 
operating range. The principle is 
the same as that behind reinforced 
concrete. 

The successful whisker-silver com- 
posite was made by causing molten 
silver to infiltrate a pack of whiskers 
in a vacuum. 

The composite, containing 12} per 
cent sapphire whiskers, held up under 
tension of 118,500 psi., which is five 
times the performance of the unrein- 
forced metal, on a strength-to-weight 
basis. 

The whiskers that went into the 
silver composite had an average 
strength of 750,000 psi. The over-all 
strength of the silver composite showed 
that the whiskers were not weakened 
by the combining process and that 
they were able effectively to bear the 
tensions. 

These whiskers were grown from 
pure aluminum using a special heating 
process. To date the experiments 
have been concentrated on combining 
the whiskers with aluminum as well 
as silver. It is planned in the near 
future to develop composites using 
steel or other high-strength refractory 
metals as a matrix. 


Computers for Evaluating Critical 
Fatigue and Alertness Levels.—Elec- 
tronic computers can be invaluable 
tools in determining critical fatigue 
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levels and alertness of astronauts, 
airmen and motorists, a Radio Cor- 
poration of America scientist told the 
Fourth International Conference on 
Medical Electronics. 

Dr. D. S. Himmelman, Administra- 
tor, Advanced Product Planning, 
RCA Electronic Data Processing Di- 
vision, said research indicates an elec- 
tronic “brain”’ can be used to excellent 
advantage in weeding out accident- 
prone individuals from among ap- 
plicants for driver’s licenses. At the 
same time, a computer can be em- 
ployed in pinpointing the critical 
fatigue levels of astronauts and com- 
mercial airline pilots. 

Dr. Himmelman said his findings 
were based on bio-medical research 
conducted in participation with Dr. 
Juan Negrin, Jr., of New York City, 
Dr. J. T. Chu of the University of 
Pennsylvania, and a group of volun- 
teer subjects. Primary purpose of 
the study was to attain evaluation, 
through the use of electronic compu- 
ters, of human capability to function 
under given stress conditions. 

Initial results of the study showed 
that in addition to space, air and 
highway safety applications, military 
personnel on vital duty—such as 
radar watch—can be monitored for 
“indications of fatigue incompatible 
with continued efficient functioning.”’ 

In adjudging the risk quotient of 
would-be motorists, Dr. Himmelman 
said, the applicant for a driver's 
license would be tested behind the 
wheel of an automobile mock-up with 
a motion picture of varying highway 
conditions replacing the customary 
windshield. Electroencephalographic 
(EEG) signals emanating from the 
subject’s brain would be transmitted 
to the electronic computer for com- 
parison with previously-determined 
standards of acceptable reaction. 


In the same manner, the RCA 
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scientist said, EEG signals from space 
vehicles could be transmitted to a 
ground-based computer system to 
analyze the alertness of spacemen. 
Commercial air line pilots could be 
tested both before and after flights. 

In the study by Dr. Himmelman 
and associates, the volunteer subjects 
were observed during periods of 
sleeplessness ranging up to 72 con- 
secutive hours, interspersed with vari- 
ous forms of visual and other stimu- 
lation. The EEG signals from the 
individual’s brain were converted to 
digital signals acceptable to an RCA 
301 electronic data processing system. 
The computer's “‘memory”’ previously 
had been stored with factors having 
a bearing on the alertness readings. 

The subjects ranged in age from 
20 to 45, and they all underwent 
the tests in a controlled environment. 

The wave forms monitored during 
the 72 hour period were related to a 
scale of scores ranging from zero to 
20 reflecting alertness or capability 
to function in humans. On the basis 
of the EEG wave forms the com- 
puter pointed to the corresponding 
alertness score of the subject under 
test. 

At present, Dr. Himmelman said, 
five points on the alertness scale are 
described by an initial prepared com- 
puter program, indicating a fatigue 
state equivalent to that induced by 
zero, 10, 20, 30, 40 and 50 hours of 
insomnia. Continued research and 
an increased sample of experimental 
subjects, will serve to improve the 
precision of the measurements, he 
said. 


Fully Automatic Camera Fits Any 


Microscope.—The first fully auto- 
matic microphotographic camera that 
fits any microscope and has both 
detail and field-integrating exposure 
mechanisms has been developed by 
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Ernst Leitz GmbH of Wetzlar, West 
Germany, and introduced to medical 
researchers and industrial laboratory 
scientists in the U. S. by E. Leitz, 
Inc., American firm which distributes 
Leitz scientific instruments and the 
Leica camera line. 

The new Orthomat micro-camera 
incorporates problem-solving features 
far more diverse in application than 
all available equipment. With a film 
magazine of 36 exposures, a researcher 
can rapidly take pictures of 36 dif- 
ferent subjects and be assured of 
perfect exposures every time. 

The extraordinary light-sensitive 
photocell is a photo multiplier tube, 
designed by R.C.A. This tube pro- 
duces a current from the available 
light, charging a capacitor to a pre- 
determined voltage. A gas discharge 
tube then starts the discharge, operat- 
ing a relay; the electrical vibration- 
free shutter opens and closes, ending 
the exposure; and the film is trans- 
ported by one frame. So automatic 
and versatile is the instrument that 
the researcher can operate the camera 
rapidly by foot pedal even as he 
watches through the microscope and 
makes adjustments with his hands. 

The Leitz Orthomat is the first 
to achieve both built-in detail as well 
as integrating measurements, oper- 
ated simply by two ‘“‘push-buttons.”’ 
An integrating measurement refers 
to an entire area. A detail measure- 
ment allows for correction of any 
particular detail within the micro- 
scopic field. 

Another outstanding feature is that 
the film magazine can be changed at 
any time, place or interval.  Ex- 
posures on color film can be changed 
at will to black and white work, and 
vice-versa, in two seconds. The self- 
contained 35-mm. camera has inter- 
changeable cassettes to hold any 
type of film desired. 
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“It has always been a problem to 
determine the exact exposure neces- 
sary for a microscopic image,” a 
company optical scientist said, ‘“‘es- 
pecially when conditions are bad, 
with dark field and tiny details. This 
instrument really takes care of 
an infinite variety of images from 
metallurgical samples to biological 
specimens.” 

The instrument is believed to be 
of profound significance to the coun- 
try’s pathologists, bacteriologists, med- 
ical technologists and others in the 
health field, as well as researchers 
in the steel, chemical, food processing, 
petroleum and a host of other indus- 
tries. Virtually eliminating the 
chance of failure because of lack of 
particular knowledge of microphotog- 
raphy or other factors is considered 
an important achievement. 

Many problem pictures are finally 
made easy. The range of exposure 
time is great. It includes those of 
photography of fluorescent specimens 
—}-hr. exposures and longer—and 
photography of living specimens, which 
have to be caught at 1/100th of a 
second. 


Fuzed Fiber Optic Plates.—Fuzed 
fiber optic plates suitable for use as 
coupling elements in image tubes have 
recently been developed by the War- 
fare Vision Branch, U. S. Army 
Engineer Research and Development 
Laboratories, Fort Belvoir, Va. 

A 1-in. diameter plate is made up 
of approximately 15,000,000 indi- 
vidual glass fibers fuzed into a solid 
piece with suitable optical insulation 
between fibers. Each fiber is only 
0.00025 in. in diameter. These plates 
are able to transfer an image from 
one side of the piece through the glass 
channels to the other side, where the 
image is reformed by a series of shaded 
dots similar to a newspaper photo- 
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graph. The resultant image is high 
in resolution and contrast with negli- 
gible distortion. A program is now 
under way to place these elements in 
image tubes and other image intensi- 
fication devices. 

This work is part of a continuing 
program in the field of micro-optics. 
Future development will include fiber 
optic image inverters, translators and 
magnifiers for use in microminiature 
image intensification devices. 


Silicon Photosensors.—Philco Cor- 
poration’s Lansdale Division has an- 
nounced a new line of silicon photo- 
sensors with high voltage capabilities 
making possible greater economies 
in the design of photo diode readout 
circuitry for data processing equip- 
ment. The new photosensors are 
capable of handling up to 20 volts 
of reverse breakdown voltage with 
under 10 ya of leakage current. They 
are designed as readout devices for 
use in conjunction with circuitry to 
convert light information into usable 
electrical signals such as pulses and 
sinusoidal waves. The photosensors’ 
high voltage, low leakage current 
capabilities represent nearly two orders 
of magnitude of improvement over 
similar devices commercially available. 

The improved characteristics offer 
a new economy to designers of data 
processing equipment using light de- 
tection as a readout method. For 
example, when used as light sensors 
for bias control in transistorized per- 
forated tape or punched card readers, 
the improved voltage capabilities of 
these diodes permit the use of direct 
coupled circuitry. Devices now in 
use, having voltage handling capa- 
bilities of little more than 2 volts, 
are not adaptable to direct coupled 
circuitry, and therefore necessitate 
the use of transformers, incurrng 
an added cost. The improvement in 
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voltage rating came about as a result 
of optimization of the diffusion process 
during production. 

The new photosensors are hermet- 
ically sealed in glass envelopes of 
0.080-in. diameter. Models are avail- 
able with active area _ positioned 
parallel to the longitudinal axis of the 
diode, to provide a side-reading capa- 
bility. Fast response time, low noise 


operation and high reliability char- 


terize the devices. A second type 
has the active area positioned per- 
pendicular to the longitudinal axis, 
and a lens for frontal exposure to 
the incident light eliminating the 
need for critical focusing. 

Though designed primarily for 
punched card and tape read-out 
systems, other applications include— 
teletypewriter coding, production flaw 
testing, signal switching, photographic 
sensing and control, optical scanning 
and sound pickup. 

Spectral response covers the range 
from visible to near infrared radia- 
tion, from 5000 angstroms to 11000 
angstroms. Peak response is at 9500 
angstroms. 


“Plastinauts” for Radiation Re- 
search.—Performance testing of three 
man-size artificial figures designed 
for biological studies of radiation 
above the atmosphere began late 
this summer by the Air Force Systems 
Command. 

The dummies are made of a rubber- 
based plastic covering a_ skeletal 
framework, and are instrumented to 
measure radiation doses in various 
parts of the body. 

These ‘“‘plastinauts’’ were  con- 
structed by Alderson Research Labo- 
ratories, Inc. of Long Island City, 
New York under a_ $35,000 con- 
tract for Systems Command's Air 
Force Special Weapons Center in 
Albuquerque, New Mexico. 
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Each dummy, which can be flown 
in any space vehicle designed to carry 
a man, weighs 162 lb. and is 5 ft. 9 in. 
tall. These dimensions are based on 
a 1950 anthropometric survey of more 
than 4,000 Air Force flying personnel. 

The dummy’s plastic covering ab- 
sorbs radiation in a manner similar 
to human tissue. Space radiation will 
be measured by small ionization 
chambers inserted in various areas 
of the figures to correspond with 
organs of the human body. Such 
measurements are important to a 
clear understanding of biological ef- 
fects of space radiation at different 
body depths. 

High-strength fabric flying suits 
have been developed for the dummies 
by the David Clark Co. of Worcester 
Mass. The suit tightly encases each 
body segment and _ stabilizes the 
dummy in the position required. 

The dummies are built to withstand 
the forces of missile launches: con- 
tinuous loadings of 10 Gs and impacts 
of 25 Gs over a wide temperature 
range from 40° F. to 140° F. 

Performance testing, conducted by 
the Biophysics Division, will include 
instrument calibration and structural 
strength studies. 

Biophysicists will consolidate data 
received from the “‘plastinaut’’ flights 
with information collected from tissue- 
equivalent dose-rate-meters, small 
pound-and-a-half instruments which 
have been carried by research rockets 
to as high as 1,600 miles. 

The dose-rate-meter uses lucite to 
simulate human tissue. 

For the dummies, a tissue-simu- 
lating plastic based on an isocyanate 
rubber, duPont Adiprene-L, is used. 
Materials like lucite would be very 
difficult to adjust to correct radio- 
logical properties in an artificial figure 
and would not have the strength, 
toughness and resilience of the adi- 
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prene-base material. The adiprene- 
base plastic has an atomic makeup 
similar to human tissue, containing 
all the trace elements found in tissue. 

Radiation-measuring instruments 
are inserted into pre-cast ducts. One 
duct is located in each arm and leg, 
one in the skull, and two each in the 
chest and abdomen areas of the 
figures. Electrical signals from the 
instruments are amplified electron- 
ically for telemetering to ground 
receiving stations. 


Automation Applied to X-Ray Crys- 
tallography.—Research into the atomic 
structure of crystals will be signifi- 
cantly accelerated by a new X-ray 
diffraction technique which is almost 
entirely automatic. Developed by 
Dr. S. C. Abrahams of Bell Telephone 
Laboratories, it involves the use of a 
high-speed, general purpose, digital 
computer to generate a control tape 
which automatically operates a special 
type X-ray diffractometer. Data 
readings from the crystallographic 
experiment are automatically recorded 
on a punched tape which is then fed 
back into the computer. The com- 
puter correlates the information, cor- 
rects for experimental factors, and 
prints out crystallographic data in a 
form convenient for analysis by a 
crystallographer. 

The new technique, called PEXRAD 
(Programmed Electronic X-Ray 
Automatic Diffractometer), was 
described in a _ paper by Dr. 
Abrahams at the American Crystal- 
lographic Association meeting in 
Boulder, Colorado. Dr. Abrahams 
said the equipment is in the final 
stages of construction and is expected 
to be in operation within a month. 

Dr. Abrahams reported that with 
PEXRAD, it will be possible to ob- 
tain more than 17,000 readings a day; 
and while they are being taken auto- 
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matically, a crystallographer will be 
left free to evaluate collected data. 
With manual methods, the most an 
experienced crystallographer can hope 
to obtain is about 3000 readings a day, 
and then probably only for one or 
two days at a stretch since the meas- 
urements require his complete 
attention. 

Dr. Abrahams said that the princi- 
ples of PEXRAD are similar to an 
automatic neutron diffractometer 
method which he and Dr. E. Prince 
developed at Bell Laboratories in 
1958. Single-crystal, automatic neu- 
tron diffractometers are currently be- 
ing used in a number of research 
laboratories which have nuclear 
reactors. 

He explained the PEXRAD method 
in detail: information such as crystal 
lattice constants, the wavelength of 
the X-ray beam, and instrumental 
constants are fed to the computer on 


punched cards, along with a compiling 


program. The computer generates a 
magnetic tape from which a punched 
paper tape is made. Fed _ into 
PEXRAD, it controls the motors 
which rotate the crystal and receiving 
counter into position, and then causes 
an X-ray beam of prescribed wave- 
length to irradiate the crystal for a 
precise length of time. The intensity 
of the beam which is diffracted by the 
crystal is measured and recorded on 
punched tape. The crystal is auto- 
matically rotated to the next position 
and the process is repeated. 

The punched tape from PEXRAD 
containing the readings is converted to 
a magnetic tape and is fed back into 
the digital computer. The computer 
integrates the intensities, corrects for 
absorption, rotation, Lorentz and 
polarization factors, and prints out 
the corrected data from which the 
location of the individual atoms can 
be calculated. 
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Dr. Abrahams said that’ with 
PEXRAD, a crystal can be oriented 
with hundredth-of-a-degree ac- 
curacy; therefore, it can be used to 
study almost any nonbiological crys- 
tal. The method could be extended 
to protein crystals, which need about 
one-thousandth-of-a-degree accuracy 
in crystal positioning. 


Memory Storage Unit for Analog 
Computers.—Electronic Associates, 
Inc., Long Branch, N. J., has an- 
nounced the first commercial sale of 
a general purpose analog computer 
equipped with a newly developed com- 
ponent that adds a high-speed memory 
storage capability. The memory com- 
ponent, called ‘‘Microstore,”’ is tran- 
sistorized and half the size of a carton 
of cigarettes. It has been under de- 
velopment and used experimentally 
for nearly a year. The commercial 
sale was made to a major chemical 
company for use in the study of 
chemical processes and associated con- 
trol systems. 

“Microstore,”” as an analog com- 
puter component, can store answers 
or values to problems for presentation 
later when the operator wants to com- 
pare new values to those previously 
computed. Two values at a time can 
be stored in one memory package for 
recall later. Analog computers, un- 
like digital, provide answers or values 
as electrical voltages that are dis- 
played graphically as curves or oscillo- 
scope displays. 

Up to 10 “‘Microstore’’ components, 
each containing two storage units, can 
be plugged into an EAI PACE 231R 
general purpose analog computer with- 
out disturbing the use of existing sum- 
ming amplifiers. According to the 
company, this design arrangement of- 
fers a distinct advantage, because the 
number of amplifiers that can be kept 
in operation determines the extent and 
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flexibility of the computer as a prob- 
lem-solving instrument. The addition 
of the memory capability to the com- 
puter permits the solution of types of 
problems that were heretofore ac- 
complished conveniently only by digi- 
tal means. 

The additional techniques made 
possible by ‘‘Microstore’’ can be ap- 
plied to a multitude of problems in all 
industries where analog computation 
applies. The most immediate and 
probably best market is in the chemi- 
cal and petrochemical fields for solving 
problems associated with various re- 
finery processes. The new component 
can be used on existing EAI large- 
scale computers. 

The company announcement also 
said that analog memory capa- 
bility was first developed by EAI for 
an Air Force computer that was de- 
livered in 1959. A second large-scale 
computer with memory capability was 
delivered earlier this year to the Na- 
tional Aeronautics and Space Admini- 
stration research facility at Langley 
Field, Va., for use in simulating vari- 
ous aspects of Project Mercury, the 
“Man in Space” program. The com- 
pany described this $1.5 million NASA 
computer as the largest general pur- 
pose analog computer ever delivered 
and said it was equipped with 45 
memory storage units. 


Electronic System for Monitoring 
Hospital Patients... omorrow’s hos- 
pital patients may have their tempera- 
ture, pulse, and breathing under con- 
stant watch from one central point on 
each hospital floor by means of a 
simple “building block’ system of 
tiny electronic measuring and radio 
transmitting devices described by 
scientists of the Radio Corporation of 
America at the Fourth International 
Conference on Medical Electronics. 
I]lustrating the concept was an 
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experimental RCA radiocardiograph 
using a sensitive transmitter only 
slightly larger than a lump of sugar. 
Placed against the skin at the finger- 
tips, on the arm, or elsewhere, the 
miniature unit detects heart action 
from the associated electrical changes 
in the body and transmits the informa- 
tion to a receiving antenna elsewhere 
in the room for display as a visible 
cardiogram on a viewing tube or an 
inked pattern on recording paper— 
thus eliminating the need for wire con- 
nections from the patient to the 
recorder. 

The application of the principle to 
routine hospital measurements was 
described by Dr. Vladimir K. Zwory- 
kin, Honorary Vice President of RCA, 
and Fred L. Hatke, of the RCA Lab- 
oratories technical staff at Princeton, 
N. J. They said that similar minia- 
turized telemetering transmitters can 
now be developed to pick up and 
transmit a variety of such measure- 
ments, including temperature and rate 
of breathing as well as heart action. 
Built into tiny packages that can be 
conveniently worn by patients, they 
said, the devices can keep a constant 
check on these routine functions and 
transmit the information to an an- 
tenna installed in a room or ward. 
From the antenna, the data would be 
carried by wires to any central loca- 
tion for display to a floor nurse or 
doctor. 

They pointed out that such meas- 
urements are now taken from time to 
time by the hospital staff as part of 
hospital routine to keep a watch on 
the general condition of each pateint. 
However, they added, “significant 
changes can occur and go completely 
undetected for some time during 
the period between _ successive 
measurements.” 

“With present methods of patient 
care, continuous monitoring of all 
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patients is impractical,” they added. 
“‘However, the simple ‘building block’ 
technique of telemetry, with parts of 
the system added in successive steps, 
can be conveniently and economi- 
cally used by all patients. And since 
all patients are in wards, semi-private, 
or private rooms, transmission is re- 
quired over relatively short distances 
from the patient to an antenna.” 

Pointing out that such a system is 
now technically feasible, the RCA 
scientists predicted that its first ap- 
plication probably will be in surgical 
recovery rooms where continuous pa- 
tient monitoring is essential. 


Quartz Clock.—A new type of clock 
which is at least ten times as accurate 
as the best conventional chronometer 
has been developed at the General 
Electric Research Laboratory. The 
clock in its present form gains or 
loses only 5 seconds a year, and even 
better performance may be obtained if 
desired. The device is small and 
rugged, and operates for about four 
months on a single size D mercury 
battery (about the size of a standard 
flashlight cell). The clock is a labora- 
tory device and no plans to produce it 
commercially have been announced. 

The new clock, which was developed 
by Robert L. Watters, operates on an 
entirely different principle from that 
of an ordinary clock. Its timing 
standard is a quartz crystal, which 
produces very stable electrical vibra- 
tions at a high frequency. The 
quartz crystal has long been recog- 
nized as an almost ideal time stand- 
ard, but its use in such a compact, 
practical chronometer became _pos- 
sible only with the recent development 
of the new solid-state electronic de- 
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vice, the tunnel diode. Previous 
clocks employing quartz crystals have 
required large amounts of electronic 
equipment to change the crystal’s 
high frequency signal to the low 
frequency required to drive a clock 
motor. They also required bulky 
storage batteries to supply the large 
amounts of power needed. 

Watters’ clock, in contrast, uses 
only three tiny tunnel diodes to con- 
vert the 100,000 cycle-per-second crys- 
tal vibrations into a 50 cycle-per-sec- 
ond current to drive the clock motor. 
This is possible because of the unique 
ability of the tunnel diode to change a 
high frequency signal into one of 
much lower frequency, something that 
is relatively difficult to do using 
vacuum tubes or transistors. In fact 
Watters developed the new clock 
originally as a method of demonstrat- 
ing the extreme stability of the tunnel 
diode in this application. Because of 
the small amount of power required 
by these electronic components, and 
a special low power clock motor pro- 
duced by General Electric's Clock and 
Timer Department, the clock can 
operate for months on a small battery. 

The tunnel diodes have proved to 
be so stable in operation that no in- 
accuracy in the clock can be ascribed 
to them. The only source of variation 
is the quartz crystal, which changes 
frequency very slightly with changes 
in temperature. In the clocks built 
so far, this inaccuracy amounts to only 
about five seconds per year over a 
wider temperature range than would 
be: expected in normal conditions of 
use. By keeping the clock at a reason- 
ably constant temperature, and per- 
haps by certain refinements in the elec- 
tronic design, even greater accuracy 
can be obtained if desired. 
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STAFF POSITIONS 


.. are available for qualified scientists and engi- 
neers in the Laboratories of The Franklin Insti- 
tute. These Laboratories comprise one of the 
country’s finest independent research organiza- 
tions, and provide R&D services to both in- 
dustry and government, in a variety of fields: 


chemistry, solid state sciences, applied physics, 
operations analysis, engineering psychology, 
and chemical, electrical, mechanical and nu- 
clear engineering. 


There are openings, in groups working on chal- 
lenging problems in both research and develop- 
ment, for men with proven abilities in any of 
these fields. 


THE FRANKLIN INSTITUTE LABORATORIES 
FOR RESEARCH AND DEVELOPMENT 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pennsylvania . .. LOcust 4-3600 


Send complete resumé to 
Mr. John E. Christ, Director of Personnel 


Or write for a brochure describing activities 
of the Laboratories 
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Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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